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An  efficient  modeling  approach  is  described  to  simulate 
the  dependence  of  remotely  sensed  radiances  at  visible  and 
near  infrared  wavelengths  (0. 4-1.1  um)  on  meteorological 
variables  such  as  relative  humidity,  wind  speed,  and  visual 
range.  Model  elements  include:  (1)  atmospheric  opacity  due 
to  Rayleigh  scattering  and  gaseous  absorption,  (2)  physical 
modeling  of  aerosol  size  distribution,  (3)  aerosol  optical 
properties  calculations,  and  (4)  implementation  of  radiative 
transfer  theory.  Sensitivity  analyses  are  performed  for  these 
variables  in  the  context  of  examining  modeled  meteorological 
influences  in  the  marine  environment  on  the  DMSP  VHR  and  LF 
sensors.  The  sensitivity  of  emergent  radiances  in  the  near 
infrared  spectral  windows  to  the  presence  of  relative  humid¬ 
ity  and  wind  speed  dependent  coarse  mode  aerosols  of  oceanic 
origin  is  demonstrated. 
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1. 


INTRODUCTION 


1 . 1  Background 

Anomalous  gray  shade  patterns  appearing  in  DMSP  (Defense 
Meteorological  Satellite  Program)  visible  VHR  (Very  High  Reso¬ 
lution)  and  LF  (Light  Fine)  imagery  are  variously  attributable 
to  low  level  haze  and  moisture,  light  fog,  cirrus  cloudiness, 
shallow  or  turbid  water,  ocean  spray,  and  sun  glint  (Fett  and 
Mitchell,  1977) .  In  a  recent  study  focusing  on  reduced  vis¬ 
ibility  situations  at  sea  (Barnes  et  al.,  1978),  DMSP  and 
concurrent  LANDSAT  MSS  (Multi-Spectral  Scanner)  data  were 
examined  to  delineate  the  roles  of  aerosol  and  water  vapor 
in  producing  anomalous  gray  shades.  Although  results  of  the 
imagery  analysis  were  not  conclusive,  an  accompanying  theo¬ 
retical  study  based  on  radiative  transfer  theory  indicated 
that  simulated  satellite-based  intensities  are  highly  depen¬ 
dent  on  the  atmosphere's  aerosol  loading  (parameterized  by 
the  surface  visual  range)  while  the  role  of  water  vapor 
(varied  over  two  orders  of  magnitude)  was  relatively  insig¬ 
nificant.  In  conducting  this  preliminary  analysis,  however, 
a  number  of  approximations  and  simplifying  assumptions  were 
adopted  including:  (a)  choice  of  a  specific  aerosol  model, 

(b)  neglecting  the  coupling  between  relative  humidity  and 
aerosol  growth,  (c)  assuming  that  the  aerosol  scattered  light 
isotropically  (independent  of  direction) ,  and  (d)  neglecting 
the  influence  of  other  meteorological  factors  such  as  wind 
speed  on  the  aerosol  size  distribution. 
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In  a  subsequent  study  relaxing  a  number  of  these  assump¬ 
tions  (Fett  and  Isaacs,  1979) ,  the  sensitivity  of  DMSP  VHR 
and  LF  sensors  in  delineating  regions  of  optically  thin 
maritime  haze  was  demonstrated  for  a  specific  case.  The 
accompanying  analysis  suggested  that  this  sensitivity  is  a 
result  of  the  system's  broad  spectral  response  extending 
into  and  peaking  at  near  infrared  (NIR)  wavelengths.  Fur¬ 
thermore,  the  calculation  presented  suggested  specifically 
that  a  number  of  factors  related  to  optical  properties  of 
maritime  hazes  including  scattering  anisotropy  and  wavelength 
dependent  extinction  are  responsible  for  the  apparent  capabil¬ 
ity  of  broadband  sensors  to  "see"  low  level  haze  while  exclu¬ 
sively  visible  sensors  cannot. 

1.2  Study  Objectives 

Since  low  level  haze  situations  with  reduced  visibility 
have  potential  adverse  impact  on  a  variety  of  Naval  electro- 
optical  systems,  these  preliminary  findings  suggested  a  re¬ 
quirement  to  further  investigate  the  capabilities  of  DMSP 
imagery  to  delineate  such  cases.  The  study  described  in 
this  document  was  undertaken  to  extend  the  investigation  of 
the  phenomena  described  above  and  to  examine  relevant  factors 
neglected  in  the  previous  work.  The  primary  objectives  of 
this  research  effort  were  to: 

(1)  Implement  a  physically  realistic  maritime  haze 
model  including  altitude  variations  of  relevant 
size  ranges  and  the  dynamic  response  of  aerosol 


size  distribution  and  complex  index  of  refraction 
to  meteorological  variables  such  as: 

.  wind  speed 
.  relative  humidity. 

(2)  Evaluate  the  wavelength  dependent  optical  proper¬ 
ties  of  the  model  maritime  haze  for  implementation 
in  the  radiative  transfer  calculations  to  follow. 
These  optical  properties  include: 

.  extinction  and  scattering  coefficients 
.  Single  scattering  albedo 
.  phase  function. 

(3)  Adopt  an  appropriate  radiative  transfer  technique 
capable  of  treating  the  mechanisms  of: 

.  multiple  scattering 

.  near  infrared  water  vapor  absorption 
.  surface  reflection 

.  scattering  anisotropy  of  the  maritime  haze 
model 

.  sun/sensor  geometry  effects. 

(4)  Compute  wavelength  and  meteorologically  dependent 
emergent  intensities  in  the  spectral  region  from 
0.4  to  1.1  um  based  on  the  previous  modeling  efforts 
and  simulate  the  response  of  the  DMSP  VHR  and  LF 
sensors  to  variations  in  meteorological  variables 

by  appropriately  weighting  results  using  the  DMSP 
sensor  response  function. 
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(5)  Investigate  the  sensitivity  of  the  DMSP  visible/ 


NIR  sensors  to  a  variety  of  representative  meteor¬ 
ological  conditions  and  subsequently  examine  the 
feasibility  of  developing  user-oriented  operational 
nomographs  based  on  these  results. 

For  the  most  part,  these  objectives  have  been  satisfied. 

1.3  Summary  of  Report 

This  document  is  organized  into  nine  sections.  Section 
2  describes  the  methodology  employed  to  accomplish  the  objec¬ 
tives  cited  above  as  well  as  an  overview  of  this  research 
effort  from  the  perspective  of  existing  optical  propagation 
models.  Sections  3  through  6  discuss  technical  implementa¬ 
tion  of  the  adopted  methodology.  Section  3  reviews  the  cal¬ 
culation  of  non-aerosol  atmospheric  optical  properties.  Sec¬ 
tion  4  discusses  modeling  the  physical  properties  of  maritime 
aerosols,  while  Section  5  includes  a  description  of  the  cal¬ 
culations  required  to  transform  from  aerosol  physical  to 
aerosol  optical  properties.  Section  6  provides  the  theoret¬ 
ical  basis  for  the  adopted  radiative  transfer  treatment. 
Results  are  presented  and  discussed  in  Section  7.  A  summary 
of  results  is  given  in  Section  8.  Finally,  Section  9  in¬ 
cludes  relevant  conclusions  and  recommendations. 
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2 .  METHODOLOGY 


2.1  Overview 

A  variety  of  mathematical  models  are  currently  avail¬ 
able  to  assess  the  influence  of  meteorological  variables  on 
in  situ,  point-to-point  atmospheric  transmission  which  may 
constrain  the  operational  effectiveness  of  proposed  or  de¬ 
ployed  electro-optical  systems.  These  codes  include  the 
atmospheric  transmittance  algorithms  developed  at  the  Air 
Force  Geophysics  Laboratories  (LOWTRAN  2,  Selby  and  McClat- 
chey,  1972;  LOWTRAN  3,  Selby  and  McClatchey,  1973;  LOWTRAN 
3B,  Selby  et  al.,  1976;  LOWTRAN  4,  Selby  et  al. ,  1978;  HITRAN, 
McClatchey  et  al.,  1973;  Rothman  and  McClatchey,  1976;  LASER, 
McClatchey  and  D'Agati,  1978),  the  optical  parameter  models 
(OPM)  adopted  by  the  Naval  Surface  Weapons  Center  (NSWC) 

(Wells  et  al.,  1977;  Hughes  and  Richter,  1979),  and  various 
battlefield  scenario  models  developed  by  the  Army  (e.g., 

Gomez  and  Duncan,  1978) .  As  described  in  the  previous  sec¬ 
tion,  this  investigation  has  focused  on  the  development  of 

computational  algorithms  to  simulate  the  parametric  response 

_2 

of  remotely-sensed,  wavelength-dependent  radiances  (mW  cm 
ym  1  sr  '*')  to  representative  variations  in  meteorological 
variables.  An  essential  step  in  the  adopted  simulation  meth¬ 
odology  is  the  evaluation  of  meteorologically-dependent  atmo¬ 
spheric  optical  properties  facilitated  by  the  application  of 
the  propagation  codes  cited  above.  This  information  is  re¬ 
quired  to  perform  radiative  transfer  calculations  to  obtain 


satellite-incident  radiances.  In  a  fundamental  sense,  there¬ 


fore,  the  models  developed  in  the  course  of  this  investiga¬ 
tion  are  an  extension  of  extant  surface  propagation  models 
and  may  provide  a  potentially  informative  link  between  avail¬ 
able  in  situ  and  remotely-sensed  data  sets. 

The  results  reported  here  are  confined  to  an  examina¬ 
tion  of  meteorological  influences  on  remotely-sensed  radi¬ 
ances  within  the  DMSP  VHR  and  LF  sensor  bandpass  illustrated 
in  Figure  1  (Fett  and  Mitchell,  1977).  This  spectral  region 
encompasses  both  visible  (0.4  to  .7  ym)  and  near  infrared 
(NIR)  (.7-1.1  ym)  wavelengths  which  arise  due  to  backscat- 
tering  of  incident  solar  radiation  from  the  atmosphere  and 
reflection  from  the  atmosphere  and  the  earth's  surface. 

The  meteorological  variables  treated  and  their  respec¬ 
tive  domains  are  given  in  Table  1.  It  should  be  noted 
that  number  density  (c)  and  visual  range  (d)  are  not  indepen¬ 
dent  variables.  In  practice,  when  visual  range  appears  as  a 
parameter,  it  has  been  used  to  compute  an  equivalent  number 
density  (see  §5.3).  Furthermore,  the  radiative  transfer 

methods  adopted  depend  on  the  product  of  number  density  and 

_2 

scale  height  [i.e.,  a  column  density  (cm  )],  not  indepen¬ 
dently  on  each.  Therefore,  increases  in  scale  height  may 
be  interpreted  as  changes  in  number  density  where  appropri¬ 
ate  (see  §5.4).  Subsequent  discussion  will  clarify  these 
points. 
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DMSP  Spectral  Bandpass 
(Fett  and  Mitchell,  1977) 
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Table  1 


Meteorological  Variables  and 
Their  Representative  Domains 


(a) 

(b) 

(c) 

(d) 

(e) 


Variable 


Relative  Humidity  (%) 
RH 

Wind  Speed  (ms 
WS 

-3 

Number  Density  (cm  ) 
N 

Visual  Range  (km) 


Aerosol  Scale 
Height  (km) 

H 


Domain 

50,  70,  80,  90,  95 
0,  5,  7,  10 
4000,  5000  (see  text 
5,  10,  23,  50 

0.5,  1.0,  1.5,  2.0 
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2 . 2  Technical  Approach 

The  general  approach  adopted  in  this  investigation  is 
summarized  in  Figure  2.  Requisite  input  parameters  include 
the  meteorological  variables:  relative  humidity  (RH) ,  wind 
speed  (WS) ,  number  density  (N)  or  visual  range  (Vr) ,  and 
aerosol  scale  heights  (H) +  and  the  geometric  configuration 
of  the  sun/atmosphere-ocean/sensor  system  as  determined  from 
the  solar  zenith  angle,  sensor  angle,  and  azimuth  angle  dif¬ 
ference.  Operationally,  these  geometric  parameters  vary 
spatially  for  individual  resolution  elements  within  a  par¬ 
ticular  sensor's  field-of-view  and  temporally  due  to  sensor 
attitude  drift.  They  are  available,  however,  from  standard 
spacecraft  ephemeris  data.  For  the  purpose  of  these  sensi¬ 
tivity  studies,  a  nadir  viewing  sensor  (sensor  zenith  angle, 
0.0°)  and  a  solar  zenith  angle  of  60°  were  assumed.  This 
provides  a  configuration  with  a  scattering  angle  of  120° 
which  is  between  those  of  the  DMSP,  NOAA5  and  SMS  spacecraft 
as  described  in  Fett  and  Isaacs  (1979).  The  results  obtained 
for  each  set  of  input  parameters  are  an  intensity  (or  radiance) 
spectrum  for  wavelengths  between  0.4  and  1.1  um,  1(A),  and  a 
DMSP  VHR  response  function  weighted  radiance,  I (DMSP). 

The  processes  linking  input  parameters  and  desired  out¬ 
put  as  illustrated  in  Figure  2  may  be  divided  into  four  basic 
modeling  efforts:  (a)  atmospheric  transmission  (i.e.,  exclus¬ 
ive  of  aerosols) ,  (b)  physical  modeling  of  aerosol  properties 

+  Depending  on  the  aerosol  model  adopted,  individual  scale 
heights  Hi  and  H2  may  be  assigned  to  individual  aerosol 
size  ranges. 
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including  relative  humidity  and  wind  speed  dependence  of 
aerosol  size  distribution  and  index  of  refraction,  (c)  de¬ 
termination  of  wavelength-dependent  aerosol  optical  proper¬ 
ties,  and  finally  (d)  radiative  transfer  theory.  To  a  great 
extent,  typical  propagation  models  (§2.1)  include  elements 
(a) ,  (b) ,  and  (c) .  Implementation  of  some  form  of  radiative 
transfer  code  is  required  to  evaluate  radiances  multiply- 
backscattered  to  a  satellite  sensor.  In  the  formalism  adop¬ 
ted,  the  radiative  transfer  theory  is  dependent  (at  each 
wavelength)  on  three  variables  which  characterize  the  ensemble 
properties  of  the  atmosphere's  optically  active  gaseous  and 
particulate  constituent  species.  These  are:  (1)  the  total 
optical  depth,  t*,  (2)  the  single  scattering  albedo,  u>0 ,  and 
(3)  the  scattering  phase  function  asymmetry  parameter,  g. 

These  three  dimensionless  variables  are  determined  from  the 
wavelength-dependent  optical  properties  of  individual  atmo¬ 
spheric  constituents  including:  (a)  scattering  by  molecules 
(i.e.,  Rayleigh  scattering),  (b)  absorption  by  molecules 
(such  as  water  vapor  and  ozone) ,  and  (c)  scattering  and  ab¬ 
sorption  by  aerosols.  Atmospheric  optical  properties  are,  in 
turn,  dependent  on  meteorological  variables  as  illustrated  in 
Figure  2.  Relative  humidity  fixes  both  the  abundance  of  water 
vapor  in  the  atmosphere  (assuming  a  model  temperature  profile) 
thereby  controlling  water  vapor  absorption  where  applicable 
(such  as  in  the  NIR  water  vapor  bands)  and  the  rate  of  aerosol 
growth  which  modifies  aerosol  size  distribution  and  prescribes 
the  appropriate  aerosol  index  of  refraction.  Wind  speed  deter¬ 
mines  the  rate  of  mechanical  addition  of  coarse  mode  (r  >  1.0  ym) 
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fugitive  sea  spray  which  in  turn  modifies  size  distribution. 
A  detailed  discussion  of  the  elements  of  the  technical 


approach  presented  above  follows  in  Sections  3  through  6. 


3.  ATMOSPHERIC  (NON- AEROSOL)  OPTICAL  PROPERTIES 

In  the  absence  of  aerosols,  atmospheric  optical  proper¬ 
ties  in  the  0. 4-1.1  ym  spectral  region  are  determined  by 
molecular  (Rayleigh)  scattering  and  absorption  by  gases  such 
as  water  vapor  and  ozone.  Both  the  normal  optical  depth  due 
to  Rayleigh  scattering,  t^,  and  that  due  to  absorption  by 
atmospheric  gases,  t  ,  have  been  evaluated  in  this  effort 
using  a  modified  version  of  computer  code  LOWTRAN  4  (Selby 
et  al.,  1978).  The  model  atmosphere  chosen  for  this  purpose 
is  the  tropical  model  given  in  Table  2. 


3.1  Rayleigh  Scattering  Optical  Depth,  r® 

The  Rayleigh  scattering  optical  depth,  t®,  is  obtained 
by  integrating  the  Rayleigh  scattering  coefficient,  B^, 
throughout  the  atmosphere: 

x®(A)  =  J&\ J(A,z)dz  (1) 


and  m  is  the  index  of  refraction  of  air  (dependent  on 

pressure,  temperature,  and  wavelength) 

•  3 

n(cm  )  is  the  local  number  density  of  air  (dependont  on 
pressure  and  temperature) 
and  A(yra)  is  the  wavelength  required. 
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Table  2 


MODEL  ATMOSPHERE  USED  AS  A  BASIS  OF  THE  COMPUTATION  OF 
ATMOSPHERIC  OPTICAL  PROPERTIES 

TROPICAL. 


lit. 

(km) 

Pressure 

(mb) 

Temp. 

<°K) 

Density 

(g/ma) 

Water  Vapor 
(g/ma> 

OtUIH' 

(g/ma) 

0 

1.  01312+03 

300.  0 

1.  167E+03 

1.  9K+01 

5.  6E-05 

1 

9.  040E+02 

294.  0 

1.  064E+03 

1.  3E+01 

5.  6E-05 

2 

8.  050E+02 

288.  0 

9.  689E+02 

9.  3E+00 

5.  4^-05 

3 

7. 150E+02 

284.  0 

8.  756E+02 

4.7E+00 

5.  IK-05 

4 

6. 330E+02 

277.  0 

7.  951E+02 

2.2E+00 

4.7E-05 

5 

5. 590E+02 

270.  0 

7.  199E+02 

1.  5E+00 

4.  5K-05 

6 

4.  920E+02 

264.  0 

6.  501E+02 

8.  5E-01 

4.  3E-05 

7 

4. 320E+02 

257.  0 

5.  855E+02 

4.7E-01 

4.  IE -05 

8 

3.780E+02 

250.0 

5.  258E+02 

2.5E-01 

3.  9E-05 

9 

3.290E+02 

244.  0 

4. 708E+02 

1.2E-01 

3.  9E-05 

10 

2.  8G0E+02 

237.  0 

4.  202E+02 

5.  0E-02 

3.  9E-05 

11 

2.  470E+02 

230.  0 

3.  74  0E+02 

1.7E-02 

4. IE-05 

12 

2. 130E+02 

224.0 

3.  316E+02 

6.  OE-03 

4.3E-05 

13 

1.  U20K+02 

217.  0 

2.  929K+02 

1.  OK -03 

4.  5K-05 

14 

1.  560E+02 

210.  0 

2.  57 8E+02 

1.  OE-03 

4.  5E-05 

15 

1.  320E+02 

204.  0 

2. 260E+02 

7.  6E-04 

4.7E-05 

1G 

1.  110E+02 

197.  0 

1.  972E+02 

6.  4E-04 

4.7E-05 

17 

9.  370E+01 

195.  0 

1.  676E+02 

5.  6E-04 

6.  9E-05 

1U 

7.  U90E+01 

199.  0 

1.  382K+02 

5.  OK -04 

9.  OK -05 

19 

6.  660E+01 

203.  0 

1.  145E+02 

4.  9E-04 

1.  4E-04 

20 

5.  G50E+01 

207.  0 

9.  515E+01 

4.  5K-04 

1.  9E-04 

21 

4.  800E+01 

211.  0 

7.  938E+01 

5.  IE -04 

2.4E-04 

22 

4.  090E+01 

215.  0 

6.  645E+01 

5.  IE -04 

2.8E-04 

23 

3.  500E+01 

217.  0 

5.  G1UK+01 

5.4K-04 

3.2E-04 

24 

3.  000E+01 

219.  0 

4. 703E+01 

6.  0E-04 

3.  4E-04 

25 

2.  570E+01 

221.  0 

4.  045E+01 

6.7E-04 

3.4E-04 

30 

1.  2201-1+01 

232.  ti 

1.  831E+01 

3.  6E-04 

2.4E-04 

35 

6.  000E+00 

243.0 

8. 600E+00 

1.  IE -04 

9.2E-05 

40 

3. 050E+00 

254.0 

4.  181E+00 

4.  3E-05 

4.  IE-05 

45 

1.  590E+00 

265.  0 

2. 097E+00 

1.  9E-05 

1.3E-05 

50 

8.  540E-01 

270.  0 

1.  101K+00 

6.  3K-06 

4.  3E-06 

70 

5.  790E-02 

219.  0 

9.  2  10E-02 

1.4E-07 

8.  GE-08 

100 

3.  000E-04 

210.  0 

5.  000K-04 

1.0E-09 

4.3E-U 

Source:  Selby  and  McCialchcy  (.1972) 
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In  practice,  the  wavelength  dependence  of  the  index  of 
refraction  may  be  incorporated  into  a  simplified  expression 
for  8^(A)  appropriate  for  standard  temperature  and  pressure 
conditions: 

e^AUkm"1]  =  1.0923xl0-3  A-4'0117  (3) 

R 

and  equation  (1)  approximated  as: 

tJ(A)  =  B|(A)  Hr  (4) 

where  HR(km)  is  the  equivalent  scale  height  for  Rayleigh  scat¬ 
tering.  For  the  tropical  atmospheric  model  in  Table  2,  this 
scale  height  equals  8.05  km. 

3.2  Absorption  Optical  Depth, 

Figure  3  illustrates  the  effect  of  absorption  by  atmo¬ 
spheric  gases  on  incident  solar  radiation  in  the  0. 4-1.1  ym 
spectral  region.  Absorbing  species  include  ozone,  O^,  oxygen, 
02,  and  water  vapor,  H20. 

The  Chappuis  bands  of  ozone  lie  between  0.45  and  0.74  ym 

in  the  visible  spectrum.  The  maximum  absorption  cross-section 

-21  2 

for  these  bands  (at  about  0.60  ym)  is  5x10  cm  yielding  an 

optical  depth  of  about  .05  (or  a  transmission  of  .95)  for  a 

19  -2 

typical  ozone  column  abundance  of  10  cm  (Goody,  1964) . 
Their  effect  on  remote  sensing  is,  therefore,  quite  small. 

The  solar  "red"  bands  of  oxygen  at  .688  and  .762  ym  are  rela¬ 
tively  strong  features,  however,  since  the  oxygen  abundance  is 


Spectral  energy  curve  of  solar  radiation  at  sea  level  and  extrapolated 
outside  the  atmosphere,  as'given  by  Pettit  (1951).  The  darkened  areas 
represent  gaseous  absorption  in  the  atmosphere. 
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constant  in  the  model  atmosphere  (as  is  ozone)  they  do  not 
vary  meteorologically.  The  most  prominent  gaseous  absorber 
within  the  spectral  region  of  interest  is  water  vapor.  Both 
visible  and  NIR  bands  appear  with  varying  strengths.  Some 
of  the  qualitative  features  of  Figure  3  can  be  delineated 
using  Table  3,  which  gives  the  quantum  mechanical  transition, 
band  center,  and  strength  of  the  water  vapor  bands  (Goody, 
1964;  McClatchey  et  al.,  1973).  (The  strength  is  defined  as 
the  integral  over  the  band  of  the  absorption  cross-section.) 
Bands  that  are  identifiable  in  Figure  3  are  marked  with  an 
asterisk.  Note  that  all  bands  are  ground-state  transitions 
(initial  state  000),  and  although  those  in  the  visible  (<  0.70 
pm)  are  relatively  weak,  the  near-infrared  bands  (p,  a,  t  ,  <|>) 
are  rather  strong. 

In  order  to  simulate  the  effect  of  variations  in  rela¬ 
tive  humidity  on  the  optical  depth  due  to  water  vapor  absorp¬ 
tion,  the  climatological  water  vapor  abundance  profile  in  the 
adopted  model  atmosphere  (Table  2)  is  replaced  in  the  lowest 
few  kilometers  by  a  constant  relative  humidity  water  vapor 
abundance  profile  given  by: 

p(Z/RH)  —  T(z)  exp  [19.85(1  t  (z)  ^ 

where  RH  is  the  desired  relative  humidity  expressed  in  percent 
and  T(z)  is  the  temperature  profile  from  Table  2.  Figure  4 
compares  the  water  vapor  abundance  profiles  for  the  climatol- 
ogically  based  LOWTRAN  4  tropical  atmosphere  and  the  constant 
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Table  3 


VISIBLE  AND  NEAR  INFRARED  WATER  VAPOR  BANDS 


'ransition 

Band  Center 

(am) 

Strength 

(an) 

000-411 

0.544 

2  x 

lO"20 

000-203 

0.572 

1  X 

10"22 

000-401 

0.592 

3  x 

10-22 

000-302 

0.592 

3  x 

l0_22 

000-321 

0.594 

2  v 

1  r>-22 

000-113 

0.632 

2  x 

IQ'23 

000-311 

0.652 

2  x 

io"22 

000-103 

0.698 

1  X 

10”21* 

000-400 

0.703 

1  X 

io-22 

000-301 

0.723 

3  x 

io-21* 

000-202 

0.723 

<2  x 

io-23 

000-221 

0.734 

6  x 

io-21* 

000-013 

0.796 

1  X 

io"22 

000-112 

0.806 

6  x 

io"21* 

000-211 

0.823 

6  x 

-23 

10 

000-210 

0.824 

-22 

1  X  10  * 

000-131 

0.847 

2  x 

10’2i* 

000-003 

0.906 

2  x 

io"21* 

000-102 

0.920 

4  x 

io'22* 

000-201 

0.942 

1  x 

io-20* 

000-300 

0.943 

6  x 

-22 

10  * 

000-121 

0.968 

2  X 

io“21* 

000-220 

0.972 

<4  x 

io'23 

000-041 

1.016 

4.8  x 

io-23 

000-012 

1.111 

1.2  x 

io-21 

000-121 

1.135 

5  x 

io-20 
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Lowtrah  4  Tropical  Atmosphere 
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WATER  VAPOR  (gm 


relative  humidity  profiles  adopted  in  this  work. 

Constant  relative  humidity  profiles  are  used  in  order  to 
maintain  consistency  between  the  specification  and  treatment 
of  optical  properties  for  water  vapor  absorption  and  those  for 
aerosol  extinction.  The  subsequently  applied  radiative 
transfer  algorothm  is  a  one  layer  model  requiring  uniform 
optical  properties.  Although,  in  a  well-mixed  boundary  layer 
water  vapor  mass  mixing  ratio  may  tend  to  be  constant  with 
height  rather  than  relative  humidity,  this  would  imply  a 
relative  humidity  gradient  given  a  representative  temperature 
profile.  Since  this  would  require  a  gradient  of  aerosol 
optical  properties,  constant  relative  humidity  profiles  are 
utilized . 

Evaluation  of  wavelength-dependent  gaseous  absorption 
optical  depths  was  accomplished  by  evaluating  transmittances , 
Tg(A,RH),  for  a  vertical  path  to  space  using  an  aerosol-free 
LOWTRAN  4  atmosphere  modified  for  constant  relative  humidity 
as  described  above  and  performing  the  operation: 

xg(X,RH)  =  -ln[Tg(A,RH) ]  .  (6) 

Although  equation  (6)  is  not  strictly  valid  for  band  models 
(i.e.,  where  the  Beer-Lambert  Law  does  not  hold),  more  soph¬ 
isticated  treatments  were  deemed  inappropriate  and  inconsis¬ 
tent  with  the  radiative  transfer  treatment  adopted. 


3.3  Total  Non- Aerosol  Optical  Depths 

The  total  non-aerosol  optical  depth,  x  ,  given  by: 

na 

T  (X,RH)  =  T  p  ( X )  +  T  (X,RH)  (7) 

na  g 

was  evaluated  for  the  wavelength  region  0. 4-1.1  ym  and  for 
relative  humidities  of  50,  70,  80,  90  and  95%  using  the  tech¬ 
niques  described  above.  As  illustrated  in  Figure  2,  the  Ray¬ 
leigh  optical  depth  (a)  and  absorbing  gas  optical  depth  (b) 
are  fundamental  intermediate  variables  in  the  subsequent 
radiative  transfer  calculation.  Figure  5  illustrates  the 
variation  of  total  non-aerosol  optical  depth  with  wavelength 
and  relative  humidity.  It  is  notable  that  at  visible  wave¬ 
lengths  (A  <  0.7  yin) ,  non-aerosol  optical  depth  is  independent 

of  relative  humidity  and  follows  the  approximate  X-4  Rayleigh 
scattering  law  given  by  equations  (3)  and  (4) .  A  barely 
perceptible  increase  in  optical  depth  in  the  vicinity  of 
0.6  ym  is  due  to  the  weak  Chappuis  band  ozone  absorption 
previously  mentioned.  At  near  infrared  wavelengths  (X  >  0.7 
ym) ,  however,  non-aerosol  optical  depth  is  dominated  by  water 
vapor  absorption  in  bands  identifiable  by  reference  to 
Table  3.  The  magnitude  of  optical  depth  near  the  centers  of 
these  bands  is  highly  dependent  on  relative  humidity.  In  the 
regions  approximately  between  bands  (0.75-0.80,  0.85-0.90, 
1.0-1.08  ym) ,  Rayleigh  scattering  is  negligible  and  non¬ 
aerosol  optical  thickness  is  near  zero. 


4.  PHYSICAL  MODELS  OF  MARITIME  AEROSOLS 


The  nature  of  the  remote  sensing  and  optical  transmis¬ 
sion  problem  in  the  atmosphere-ocean  system  is  extensively 
modified  by  the  presence  of  marine  aerosols.  These  perva¬ 
sive  particulates  represent  a  challenging  problem  to  radia¬ 
tive  transfer  modelers,  due  to  their  inherent  spatial  and 
temporal  variability  and  the  apparent  dependence  of  their 
physical  properties  on  meteorological  variables  such  as 
relative  humidity  and  near-surface  wind  speed.  Reference  to 
Figure  2,  indicates  that  aerosol  optical  properties  (deter¬ 
mined  by  Mie  theory  calculations)  are  dependent  on  the 
physical  properties  of  the  aerosol  including  size  distribu¬ 
tion  and  refractive  index. 

In  this  section,  the  physical  mechanisms  determining 
the  relationship  between  meteorological  variables  and 
aerosol  physical  properties  are  elucidated  and  the  specific 
physical  models  adopted  in  this  work  are  described. 

4.1  Aerosol  Size  Distribution 

The  fundamental  link  between  meteorological  variables  and 
marine  aerosol  optical  properties  is  the  marine  aerosol  size 
distribution,  n(r) [units  cm“^ym”^].  The  size  distribution 
expresses  the  number  density  of  aerosol  "droplets"  of 
radius  r  within  a  size  range  from  r  to  (r+dr) .  Thus,  the 
total  number  density  of  aerosol  droplets  in  all  size  ranges 
is : 

00 

N  (cm-3)  =  y  n(r)  dr  (8) 

o 
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Extensive  measurements  have  indicated  that  the  aerosol 
number  density  distribution  with  size  range,  n(r)  can  be 
reproduced  by  a  mathematical  model  which  assumes  a  finite 
number  of  modes  or  size  ranges  (Whitby  and  Sverdrup,  1978) , 
generally  three , consisting  of,  in  order  of  increasing 
radius:  (a)  Aitkin  nuclei,  (b)  an  accumulation  mode  (mass 

median  radius  0.1  to  1.0  ym) ,  and  (c)  a  coarse  mode  (mass 
median  radius  greater  than  1.0  ym)  .  Since  the  effect  of 
Aitkin  nuclei  on  optical  properties  in  the  0.4  to  1.1  ym 
range  is  negligible  (Nilsson,  1979),  they  are  not  treated 
in  subsequent  discussions.  It  is  generally  convenient  to 
describe  each  of  the  two  remaining  modes  by  a  tractable 
analytical  expression.  This  is  particularly  advantageous 
when  size-fractionated,  compositional  data  suggest  [as  in 
the  case  of  maritime  aerosol  over  remote  oceanic  regions 
(Meszaros  and  Vissy,  1974) ]  that  the  accumulation  mode  and 
coarse  mode  have  distinct  source  mechanisms  and  compositions 
Although  many  functional  forms  have  been  assumed  in  the 
past,  including  Junge  (1963;  1972)  and  Deirmendjian  (1969) 
distributions,  a  particularly  convenient  representation  is 
the  log-normal  size  distribution  (Dennis,  1976): 

n(r)  =  N  [exp (-ln2r/rn/2  ln2a)  ]/r/2it  Ino  (9) 

with  number  median  radius  r  and  standard  deviation  a.  The 

n 

total  number  density  of  particles  is  normalized  to  N  (cm_3) . 
On  a  leg  probability  plot,  the  number  median  is  that  value 
with  50%  of  the  distribution's  number  density  above  and  below 
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The  distribution  is  often  given  in  terms  of  the  mass  or 
volume  median  radius  rm  related  to  the  number  median  radius 
through  the  standard  deviation  by: 

In  2rm  =  In  2rn  +  3.0  ln^cr.  (10) 

The  maritime  aerosol  optical  properties  model  adopted  by  the 
Air  Force  (Shettle  and  Fenn,  1976;  1979)  and  implemented 
within  their  standard  atmospheric  transmission  codes  (Selby 
et  al.,  1978)  are  formulated  based  on  the  number  mean  radii 
and  standard  deviations  [with  reference  to  equation  (9)  above] 
shown  in  Table  4.  One  advantage  of  this  bimodal  log-normal  pre¬ 
sentation  is  its  applicability  to  incorporation  within  optical 
parameter  evaluation  models  which  include  the  effects  of 
relative  humidity  on  aerosol  growth  (Nilsson,  1979;  Shettle 
and  Fenn,  1979) .  The  gross  effect  of  a  positive  (negative) 
change  in  relative  humidity  on  an  aerosol  size  distribution 
initially  in  equilibrium  with  its  environment  is  an  increase 
(decrease)  in  the  average  size  of  the  aerosol  droplets 
within  the  distribution.  The  magnitude  of  the  size  change 
is  dependent  on  the  initial  size  and  composition  of  the 
aerosol.  The  bimodal  log  normal  representation  permits  the 
number  median  radius  for  each  mode  to  be  suitably  modified 
for  aerosol  growth  while  keeping  the  total  number  of  parti¬ 
cles  (mode  number  density)  constant. 

This  allows  for  investigation  of  variation  in  optical 
properties  with  changes  in  relative  humidity  assuming  that 
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aerosol  production  and  loss  mechanisms  remain  constant. 

For  example.  Table  5  presents  relative  humidity  dependent 
number  median  radii  for  the  bimodal  maritime  model  described 
in  Table  4.  The  specific  rate  of  growth  with  changes  relative 
humidity  depends  on  the  aerosol  growth  law  adopted  (see 
Section  4.3)  for  each  mode. 

For  example.  Table  5  suggests  that  the  increase  in 
median  size  between  50  and  95%  relative  humidities  is  a 
factor  of  about  1.75  greater  in  the  coarse  mode  than  in  the 
accumulation  mode.  The  relative  humidity  dependent  size 
distributions  given  by  the  parameters  in  Table  5  using  the 
log  normal  size  distribution  [equation  (9) ]  are  illustrated 
in  Figure  6.  (accumulation  mode+)  and  Figure  7  (maritime 
model-both  modes) ,  respectively,  assuming  constant  standard 
deviations  (as  given  in  Table  4)  for  each  mode.  One  disad¬ 
vantage  of  the  siz.e  distribution  models  cited  above  with 
respect  to  the  present  study,  is  the  neglect  of  the  dynamic 
response  of  the  sea  spray  produced  component  of  the  maritime 
aerosol  (i.e.  the  coarse  mode)  to  variations  in  prevailing 
wind  speed.  Measurements  suggest  that  the  number  density 
of  sea  spray  produced  aerosol  droplets  increases  with 
increased  wind  speed,  especially  the  larger  particles 
(Junge,  1960) . 

The  maritime  model  accumulation  mode  properties  are  equiva¬ 
lent  to  those  of  the  tropospheric  model  in  Shettle  and  Fenn 

1979. 
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Size  Distribution 
Parameters 
(log  normal) 


Aerosol  Model  N(cra"3)  rn(ym)  o 


Type 


Maritime 

Rural 
Aerosol 
Mixture 
(water  soluble 
Aerosols) 


Continental 

Origin  1.0  .03  2.24 

(accumulation 
mode) 


Marine 

Origin 

(coarse 

mode) 


1.0  0.3  2.51 


Sea  Salt 
Solution 


Table  4.  Parameters  of  maritime  aerosol  size  distributions 

used  in  AFGL  LOWTRAN  models  (Shettle  and  Fenn,  1979) . 


Number  Mode  Radius 

(ym) 

Relative 

Humidity 

Accumulation 
%  (o  =  2.24) 

Coarse 
(o  =  2.51) 

50 

.02748 

.1711 

70 

.02846 

.2041 

80 

.03274 

.3180 

90 

.03884 

.3853 

95 

.04238 

.4606 

Table  5. 

Humidity  Dependent  Mode 
(Shettle  and  Fenn,  1979) 

Radii  Maritime  Model 

Size  Distribution  (cm 


■1 


Figure  6 

Maritime  Model  (accumulation  mode  only) 
[from  Shcttle  and  Fonn,  1979] 


Figure  7 


Maritime  Model  (accumulation  and  coarse  modes) 
(from  Shettle  and  Fenn,  1979] 


An  alternative  size  distribution  model  incorporating 
this  dependence  adopted  by  the  Naval  Surface  Weapons  Center 
(NSWC)  utilizes  a  linear  combination  of  Junge  and  Deirmendjian 
(or  modified  D  distributions  for  the  accumulation  and 
coarse  modes,  respectively.  Based,  in  part,  on  an  earlier 
model  by  Wells  et  al.  (1977)  and  modified  by  Katz,  the  size 
distribution  for  the  Munn-Katz  model  is  given  by: 


n(r) 


1.7(|)"4+  1.62(C1+C2V6)exp[^  -  8 . 5  (|)  r  )F~1  (|)  , 

(11) 


where  n(r) 


=  number  of  particles  per  cm 


3 


per  ym 


(radius) 


radius  (ym) 
altitude  (m) 

scale  height  set  at  800  m  for  Z  <  1  km 


0.81  exp 


r  .  066S  , 
l1.058-SJ 


S  = 
V  = 


F  = 

r  = 

when  V  _<  7 
V  >  7 


saturation  ratio  (relative  humidity  *  100) 
wind  factor  scaled  with  surface  wind  VQ 
.  5  m/s  for  0  <_  VQ  <_  4  m/s 
(VQ  -  3.5)  m/s  for  VQ  >  4  m/s 

1  +  <Fo>3 

0.384  -  .00293V1,25 

m/s,  C1  =  350,  C2  =  103,  6  -  1.15 

m/s,  Cx  =  0,  C2  =  6900,  6  =  0.29 


In  this  model,  the  coarse  mode  aerosol  [second  term  in  (11)] 
is  dependent  on  prevailing  wind  speed,  due  to  contributions 
of  fugitive  sea  spray  components  to  this  size  range.  This 
dependence  follows  the  data  of  Woodcock  (1953)  .  Figures  8 
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SIZE  DISTRIBUTION  (cm "Vm"1) 


RADIUS 


Figure  8 

Dependence  of  NSWC  model  size  distribution 
on  wind  speed  (z  =  Om,  RH  =  80%) 
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Figure  9 

Dependence  of  NSWC  model  size  distribution 
on  relative  humidity  (z  =  Om, WS  =  Oms-1) 
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and  9  illustrate  the  dependence  of  the  size  distribution 
given  by  Equation  (11)  on  wind  speed  and  relative  humidity, 
respectively.  The  relative  humidity  dependence  given  by  the 
growth  factor  F  in  (11)  used  in  both  modes  is  that  due  to 
Fitzgerald  (1975,  1979). 

A  recent  comparison  between  LOWTRAN  (Selby  et  al., 

1978)  and  NSWC  codes  has  been  performed  (Hughes  et  al . , 

1979)  ,  which  illustrates  their  validity  for  predicting 
visible  and  near  infrared  extinction.  It  may  be  argued  that 
due  to  compositional  and  size  distinctions  between  the  accum¬ 
ulation  and  coarse  modes,  different  growth  factors  should  be 
used.  An  advantage  of  the  coarse  mode  expression  in 
equation  (11)  above  is  that  the  number  density  realistically 
increases  with  wind  speed  as  new  particles  are  mechanically 
added.  This  is  illustrated  in  Table  6. 

4 . 2  Effect  of  Relative  Humidity 

Relative  humidity  directly  affects  both  the  equilibrium 
size  (i.e.  radius)  of  an  aerosol  droplet  and  the  complex 
index  of  refraction.  One  approach  to  including  relative 
humidity  (RH)  effects  is  based  on  evaluating  the  humidity 
growth  factor  F(f)  defined  as: 

F  ( f )  =  (12) 

ro 

where  r(f)  is  the  equilibrium  radius  of  the  particle  at 
f  =  (RH/100%)  and  rQ  is  the  radius  of  the  dry  particle. 

Given  the  dry  particle  size  distribution  n(rQ),  the  response 
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to  relative  humidity  is  obtained  by  substituting  the 
appropriate  growth  factor,  such  as: 

n [r (f ) ]  =  n[F(f)r0]  (13) 

where  n[r(f)]  is  the  size  distribution  at  RH=f .  This  assumes 
that  the  total  number  density  is  constant.  Alternatively, 
if  the  number  density  distribution  is  specified  by  a  standard 
size  distribution  function  (such  as  log  normal,  for  example), 
the  relative  humidity  effects  can  be  obtained  by  multiplying 
the  appropriate  median  radius  parameter  by  the  growth  factor. 

The  complex  refractive  index  at  relative  humidity  f, 
mf/  is  obtained  by  volume  mixing  the  refractive  index  of  dry 
particles  mQ  and  that  of  water  ni,  according  to: 

mf  =  moF(f)~3  +  mw[l  -  F(f)'3]  (14) 

This  combined  approach  has  been  used  in  most  modeling  of 
relative  humidity  effects  on  aerosol  size  distribution 
(Nilsson,  1979;  Shettle  and  Fenn,  1979;  Hughes  and  Richter, 
1979)  . 

The  appropriate  growth  factor  r/rQ=F  is  evaluated  based 
on  aerosol  growth  theory.  The  theoretical  basis  of  aerosol 
growth  has  been  discussed  by  numerous  authors  (Nilsson,  1979; 
Kohler,  1926;  Hanel,  1968,  1972;  Junge  and  McLaren,  1971; 
Winkler,  1973;  Fitzgerald,  1975,  1979,  among  others).  In 
general,  a  particle  model  is  adopted.  The  average  aerosol  is 
assumed  to  be  an  internal  mixture  (Winkler,  1973)  made  up  of 


several  specifiable  substances  mass-weighted  according  to 
their  total  proportion  in  the  distribution.  At  zero  relative 
humidity,  the  particles  are  assumed  to  consist  of  water- 
soluble  fraction  mass  E  and  a  complementary  insoluble  mass 
fraction,  (1  -  E) .  As  relative  humidity  increases  a  fraction 
E,  the  water-soluble  part  gradually  dissolves,  and  the  aero¬ 
sol  droplet  grows  by  absorption  of  water  vapor  (Fitzgerald, 
1975) .  The  mass  of  the  dry  particle  is  given  by: 

mo  =  aPoro3  (15> 


where  rQ  is  the  radius  of  the  dry  particle,  pQ  is  the  mean 
density  of  the  dry  particle,  and  a=4n/3.  The  mean  density 
of  the  dry  particle  p0  is  defined  by  summing  soluble  and 
insoluble  volumes  to  give: 


Pn  — 


pspi 


Ep.  +  (1-E)  p 
i  s 


(16) 


where  pg  is  the  density  of  the  water  soluble  fraction,  and 
p.^  is  the  density  of  the  insoluble  fraction.  The  mass  of 
the  water  soluble  fraction  mdg  is: 


mds  -  Emo 


(17) 


while  the  mass  of  dissolved  matter  is: 


mg  =  e (f)mds  =  e (f)EmQ 


(18) 


where  the  dissolved  fraction  e  depends  on  relative  humidity. 
The  mass  of  water  in  a  particle  of  size  r  in  equilibrium 


with  RH=f  is: 


mw  =  a(r3-r03)pw 

=  a[F(f) 3-l]r03  pw 


(19) 


where  pw  is  the  density  of  water.  Based  on  (15)  -  (19), 
the  density  of  the  salt  solution  p '  and  the  molecular  weight 
of  the  salt  solution,  M'  may  be  derived  by  volume-weighting: 


P '  = 


psg(r) 


m 


ms  +  mw  = 
s/ps+  mv/pw  r3pg+  r3(eEpQ-ps) 


M' 


ms  +  mw 
™M  +  m  /M 


MsMw9  (r) 


V  l*W  M  ^ m  c'lPn  — m  \  v  3 


(mweEp0-ms)r^ 


where : 


g  (r)  =  r3+(£EpQ-l)r3  . 


(20) 


(21) 


(22) 


and  Mg,  are  the  molecular  weights  of  solute  and  water, 
respectively. 

The  growth  factor  may  be  evaluated  from  the  relationship 
between  relative  humidity  and  the  equilibrium  radius  of  an 
aqueous  solution.  The  ratio  of  the  equilibrium  vapor 
pressure  over  the  surface  of  the  aerosol  droplet  of  radius 
r,  P  ,  and  that  over  a  planar  water  surface,  Pm  is  given  by: 

L  oo 


exp 


r2o'M' 

lrp'RvT 


I  Cl  + 


‘Vs,-! 


Mm 
s  w 


(23) 
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where  a '  = 


aw<T)  = 
b  = 

M'  = 

r  = 

p'  = 

Rv  = 
T  = 

i  = 


M  = 
w 

M  = 
s 


m  = 
w 


m  = 
s 


surface  tension  of  the  solution 
ow(T)  +  b(ms/mw) 

surface  tension  of  water  at  temperature  T 

a  constant  depending  on  solute 

molecular  weight  of  the  solution 

radius  of  the  droplet 

density  of  the  solution 

specific  gas  constant  of  water  vapor 

temperature  of  the  droplet  (K) 

van't  Hoff  factor,  which  is  determined  empiri¬ 
cally,  and  indicates  the  degree  of  dissociation 

molecular  weight  of  water 
molecular  weight  of  the  salt 
mass  of  water,  and 
mass  of  salt 


Assuming  instantaneous  equilibrium  with  surrounding  air 
(i.e.  time-dependent  case) ,  the  ratio  above  (23)  is  identi¬ 
cally  equal  to  the  relative  humidity,  or  substituting  from 
(15)  -  (22)  : 


f  =  exp 


a  (T)+b 

(r/rc 


eEPp  ' 
Vr0)3-1 


MsMw[(r/ro)3~1+£Epo] 
Ms[ (r/r0)3-l)  +MweEpoJ 


r0(r/r0) 


[ (r/rQ) 3-l+eEpQ]  ps 
°  [ (r/r0) 3-l]p  +eEp0 


-RvT 


(24) 


+  -Mw  eEpo 
1  +  i  n- 


Ms  (r/rn) 3-l 


-1 
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where  the  growth  factor  F(f)  is  r/rQ  in  the  expression  above. 
Figure  10  (Fitzgerald,  1979)  illustrates  a  comparison 
between  a  number  of  theoretical  models  of  F(f)  and  measure¬ 
ments  for  maritime  aerosols.  In  this  effort  we  have  sought 
to  distinguish  between  growth  factors  for  the  accumulation 
(or  fine)  mode  and  for  the  coarse  mode.  Several  candidate 
calculations  for  each  of  these  modes  are  presented  in 
Figures  11  and  12,  respectively.  Note  that  for  a  given 
relative  humidity  in  the  range  50-95%,  the  growth  factor  for 
the  coarse  mode  is  generally  greater  than  that  for  the  fine 
(accumulation)  mode.  Above  80%  relative  humidity,  the 
coarse  mode  growth  factors  of  Fitzgerald  (1975,  1979)  and 
Shettle  and  Fenn  (1979)  on  approximately  the  same.  Below 
80%,  however,  Fitzgerald's  values  essentially  average  over 
the  hysteresis  phenomenon  between  40  and  76%  relative 
humidity  while  the  others  inherently  assume  the  onset  of 
deliquessence  at  about  70%. 


4.3  Candidate  Physical  Models 

Based  on  the  previous  discussion,  two  models  were  adopted 
to  perform  the  sensitivity  studies  relating  meteorological 
variables  to  emergent  intensities  required  in  this  effort. 

The  first  of  these  is  the  Shettle  and  Fenn  (1979)  relative 
humidity  dependent  maritime  model  prescribed  by  the  param¬ 
eters  in  Table  5  and  the  log  normal  distribution  for  each 
mode  given  by  equation  (9). 
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Figure  11 

Growth  Factor  (accumulation  mode) 
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Figure  12 

Growth  Factor  (coarse  mode) 
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The  second  model  is  a  hybrid  expression  consisting  of  a 
linear  combination  of  the  accumulation  or  tropospheric  mode 
from  Shettle  and  Fenn  (1979)  and  the  wind  speed  dependent 
coarse  mode  from  the  Munn-Katz  model.  In  employing  this 
model,  the  aerosol  growth  factor  [equation  (12)]  for  the 
accumulation  (fine)  mode  is  that  used  by  Shettle  and  Fenn 
(1979)  while  that  in  the  coarse  mode  is  that  given  by  Fitz¬ 
gerald  (1975,  1979)  [see  Figures  11  and  12].  Two  advantages 
of  this  approach,  therefore,  are  that:  (a)  aerosol  growth 

rates  with  relative  humidity  are  distinct  between  the  modes 
and  (b)  wind  speed  dependence  is  treated  in  the  coarse  mode. 

In  employing  this  hybrid  model,  number  densities  in  the  coarse 
mode  vary  with  relative  humidity  and  wind  speed  as  given  in 
Table  6,  while  the  accumulation  mode  number  density  is  either 
held  constant  or  normalized  to  correspond  to  a  specified  sur¬ 
face  visual  range.  Figures  13 (a) -(d)  illustrate  the  hybrid 
size  distribution  variation  with  relative  humidity  at  fixed 

wind  speeds.  The  number  density  in  the  fine  mode  is  set  at 
_  3 

5000  cm  .  Figures  14(a)  and  (b)  illustrate  the  variation 
of  the  same  distribution  with  wind  speed  (0,  5,  7,  10  ms-1) 
for  fixed  relative  humidities  of  50  and  95%,  respectively. 
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Figure  13(a) 

Hybrid  Size  Distribution  Variation 
with  Relative  Humidity  (WS  =  0.0  ms-1) 
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Figure  13(b) 

Hybrid  Size  Distribution  Variation 
with  Relative  Humidity  (WS  =  5  ms'1) 
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Figure  13(c) 

Hybrid  Size  Distribution  Variation 
with  Relative  Humidity  (WS  =  7  ms“^-) 
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Figure  13(d) 

Hybrid  Size  Distribution  Variation 
with  Relative  Humidity  (WS  =  10  nic'M 
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Figure  14 (a) 

Hybrid  Size  Distribution  Variation 
with  Wind  Speed  (RH  =  50%) 
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Figure  14 (b) 

Hybrid  Size  Distribution  Variation 
with  Wind  Speed  (RH  =  95%) 
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5.  AEROSOL  OPTICAL  PROPERTIES 

5.1  Mie  Theory  Calculations 

Given  the  dependence  of  size  distribution  on  wind  speed 
and  relative  humidity  (i.e.  growth  factor)  and  model  aerosol 
complex  index  of  refraction  on  growth  factor  [see  equation 
(14)]  classical  Mie  theory  (Mie,  1908)  may  be  used  to  eval¬ 
uate  the  optical  properties  required  for  radiative  transfer 
calculations  (see  Figure  2) .  Details  of  the  Mie  theory  may 
be  found  in  a  number  of  standard  references  Deirmendjian, 
1969;  van  de  Hulst,  1957;  Hansen  and  Travis,  1974;  McCartney, 
1976)  . 

Specifically,  these  optical  properties  include:  (a) 
extinction  coefficient,  Be;  (b)  single  scattering  albedo, 
a)0;  and  (c)  angular  scattering  function  P(0)  or  asymmetry 
factor,  g.  Each  of  these  parameters  is  a  function  of  funda¬ 
mental  variables  such  as  wavelength.  A,  size  distribution, 
n(r),  and  complex  index  of  refraction,  in.  The  extinction 
coefficient,  6  (A),  for  an  aerosol  with  size  distribution, 
n(r),  is  given  by: 

( A)  =  Tr/r2Qe(r,A,m)n(r)dr  (25) 

o 

where  Qe  is  the  Mie  efficiency  factor  for  extinction  for  a 
particle  of  radius  r.  Similar  expressions  are  applicable 
for  the  absorption  coefficient,  B  (A),  and  scattering 

a 
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coefficient,  B  (A),  substituting  Mie  efficiency  factors  Q 
and  Qs ,  respectively  (where  Qe  =  Qa  +  Q  ) .  These  efficiency 
factors  are  computed  using  readily  available  computer  pro¬ 
grams.  (Hansen  and  Travis,  1974;  Dave,  1969;  1970). 

The  single  scattering  albedo  is  the  ratio  of  scattering 
coefficient  to  extinction  coefficient  for  a  given  aerosol 
size  distribution.  Based  on  the  discussion  above: 

w0(A)  =  (3  U)/3  (A)  =  JrQ  n(r)//r  Qen(r)dr.  (26) 

The  angular  scattering  (or  phase)  function  determines 
the  fraction  of  incident  light  at  a  particular  wavelength  A 
scattered  from  the  direction  of  the  source  into  the  direc¬ 
tion  0  degrees  (0<G<180)  from  the  source.  The  expression 
for  the  angular  scattering  function  is: 

2  °° 

P  (0)  =2W~/[il(0)  +  i2(e)  ]n(r)dr  (27) 

where  i^  and  ^  are  the  Mie  angular  intensity  functions  for 
a  particulate  of  radius  r  and  complex  of  refraction  m  at  a 
specific  wavelength  A.  These  functions  can  be  evaluated 
using  available  Mie  theory  computer  codes.  For  a  number  of 
approximate  radiative  transfer  techniques  such  as  the  one 
described  in  Section  6,  it  is  useful  to  characterize  the 
aerosol  angular  scattering  function  by  a  single  asymmetry 
parameter,  g,  defined  as  (Joseph  et  al,  1976): 
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g 


cosG  P(Q)  d(cos0) 


(28) 


Values  of  the  g  factor  range  from  +1.0  for  complete  forward 
scatter  to  -1.0  for  complete  backscatter. 

5.2  Meteorological  Dependence  of  Aerosol  Optical  Properties 
Mie  theory  calculations  of  Be(X),  w0(X),  and  g(X)  were 
required  for  the  model  aerosol  size  distributions  described 
in  section  4.3  within  the  wavelength  region  0. 4-1.1  ym  and 
for  representative  values  of  relative  humidity  and  wind  speed. 
Necessary  values  for  the  maritime  model  and  accumulation  mode 
of  the  hybrid  aerosol  model  were  available  from  tabulations 
in  Shettle  and  Fenn  (1979).  For  the  coarse  mode  of  the 
hybrid  model,  Mie  theory  calculations  were  performed  using  a 
modified  version  of  the  code  from  Hansen  and  Travis  (1974) . 
Appropriate  wavelength  dependent  indices  of  refraction  were 
obtained  by  volume  weighting  according  to  equation  (14) using 
Fitzgerald's  (1975,  1979)  growth  factor.  Indices  of  refrac¬ 
tion  for  water  and  sea  salt  were  obtained  from  Hale  and 
Qerry  (1973)  and  Volz  (1972),  respectively,  as  presented  in 
Shettle  and  Fenn  (1979) . 

5.2.1  Maritime  Model 

Extinction  coefficients  (km  ^)  and  asymmetry  parameters 

_  3 

for  the  maritime  model  (assuminq  N=4000  cm  )  as  a  function 
of  relative  humidity  in  the  spectral  range  0. 4-1.1  ym  are 
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presented  in  Figures  15(a)  and  (b) ,  respectively.  Note 
that  both  extinction  coefficient  and  asymmetry  parameter  in 
this  model  are  relatively  wavelength  independent  within  this 
spectral  interval,  but  strongly  influenced  by  relative 
humidity. 

5.2.2  Hybrid  Model 

Figures  16(a)  and  (b)  illustrate  the  wavelength/relative 
humidity  dependence  of  the  extinction  coefficient  (km  ^)  and 
asymmetry  parameter,  respectively,  of  the  hybrid  model  fine 
mode  only  (N=5000  cm  ^) .  For  this  mode,  extinction  decreases 
with  increasing  wavelength  with  a  dependence  which  may  be 
approximated  by: 

6  (A)  «  const  /\a  (29) 

where  the  Angstrom  (1961)  coefficient  is  such  that 
1.0  _<  a  <_  2.0.  Additionally,  as  wavelength  increases 
scattering  asymmetry  decreases. 

Analogous  calculations  for  the  hybrid  model,  wind  speed 
dependent  coarse  mode  are  presented  in  Figures  17 (a, b), 
18(a,b),  19(a,b)  and  20(a,b)  for  wind  speeds  of  0,5,7,  and 
10ms  ^ ,  respectively.  Generally,  although  the  magnitude  of 
the  extinction  due  to  the  coarse  mode  is  less  than  that  in 
the  fine  mode  [Figure  16 (a) J,  the  coarse  mode  extinction  is 
much  less  wavelength  dependent  and  much  more  dependent  on 
relative  humidity  due  to  its  larger  sized  particles.  For 
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Hybrid  Model  Fine  Mode 
(from  Shettle  and  Fenn,  1979) 
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Asymmetry  Parameter-Hybrid  Model  Fine  Mode 
(from  Shettle  and  Fenn,  1979) 
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Asymmetry  Parameter-Hybrid  Model 
Coarse  Mode  (WS  *  0  ms-*) 
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Figure  18(b) 

Asymmetry  Parameter-Hybrid  Model 
Coarse  Mode  (WS  =  5  ms-1) 
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Figure  19(a) 
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Figure  19(b) 

Asymmetry  Parameter-Hybrid  Model 
Coarse  Mode  (WS  =  7  ms"-1-) 
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the  same  reason,  its  range  of  asymmetry  parameters  is 
greater  than  that  in  the  fine  mode. 

5.3  Number  Density/Visual  Range  Normalization 

As  remarked  earlier,  the  overall  normalization  of  the 
aerosol  size  distributions  adopted  may  be  discussed  either 
in  terms  of  number  density  (cm  3)  or  visual  range  (km) . 
Aerosol  extinction  at  0.55  ym,  Be  (0 . 55),  and  visual  range, 

VR,  are  related  approximately  by  the  relation  (Koschmieder , 
1924;  Horvath  and  Noll,  1969): 

VR(km)  =  3 . 912/ [ ge (0 . 55)  +  gj(0.55)]  (30) 


where  the  Rayleigh  scattering  coefficient  at  0.55  ym  given 
by  equation  (3)  is  0.012  km  . 

Using  equation  (30)  and  the  aerosol  extinction  coeffi¬ 
cients  evaluated  in  §5.2,  corresponding  visual  ranges  may  be 

•  -3 

calculated  assuming  the  fixed  number  densities  of  4000  cm 

_  3 

for  the  maritime  model  and  5000  cm  for  the  fine  mode  of  the 
hybrid  model.  For  the  hybrid  model  including  the  Munn-Katz 
coarse  mode,  the  coarse  mode  number  densities  are  not  con¬ 
stant  but  are  given  in  Table  6.  The  visual  ranges  corres¬ 
ponding  to  the  adopted  models  are  presented  in  Table  7  as  a 
function  of  relative  humidity  and,  where  appropriate,  wind 
speed . 

Alternatively,  it  may  be  desirable  to  consider  constant 
visual  range  scenarios.  In  this  case,  it  is  necessary  to 
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adjust  the  total  number  density  in  the  distribution  to 
provide  the  appropriate  value  of  aerosol  extinction  at  0.55 
Vim.  Number  densities  adjusted  to  correspond  to  visual 
ranges  of  5,10,23,  and  50km  are  presented  in  Table  8.  For 
the  hybrid  model  it  was  decided  to  retain  the  coarse  mode 
number  densities  given  in  Table  6  and  adjust  only  the  fine 
mode  number  densities. 

For  the  most  part  subsequent  simulated  radiance  results 
are  presented  in  terms  of  the  fixed  number  density  models 
as  given  in  Table  7.  However,  a  few  results  are  presented 
as  functions  of  visual  range. 

5.4  Aerosol  Optical  Depths,  xe,s 

i 

Aerosol  optical  depths  for  scattering  and  extinction 
required  in  the  radiative  transfer  calculation  may  be 
evaluated  based  on  the  scattering  and  extinction  coefficients 
described  above  for  each  adopted  model.  The  scattering  or 
extinction  optical  depth  for  a  given  mode  i  is  defined  in 
analogy  to  equation  (1)  as: 

t l'S  (X)  =  S$l's  ( X, z )  dz  (31) 

If  an  approximately  exponential  dependence  of  aerosol  number 
density  with  altitude  is  assumed  with  scale  height  (or  a 
uniformly  mixed  situation  with  mixing  height  H^) ,  equation 
(31)  may  be  approximated  as  : 


(32) 


r 


T®'S  (X)  =  I®,S  (X)  H± 

_0  g 

where  3^'  represents  a  near  surface  or  uniformly  mixed 
value.  In  general,  subsequent  results  assume  a  scale  height 
of  1.2  km  for  the  maritime  model  and  fine  mode  of  the  hybrid 
model  based  in  part  on  Elterman  (1970)  .  For  the  Munn-Katz 
coarse  mode  a  scale  height  of  0.8  km  is  used  (Hughes  and 
Richter,  1979) .  These  choices  of  scale  height  qualitatively 
describe  the  relatively  longer  lifetime  of  accumulation  mode 
aerosols  due  to  their  smaller  sizes  and  resultant  resistance 
to  loss  due  to  gravitational  settling.  For  the  sensitivity 
studies  which  are  performed  on  scale  height,  however,  it  should 
be  noted  that  the  scale  height  value  may  be  considered  as  the 
product  of  a  physical  height  and  a  dimensionless  constant 
which  increases  or  decreases  the  total  number  density  (with 
associated  changes  in  surface  visual  range).  Therefore,  these 
calculations  may  be  viewed  as  constant  scale  height  calcula¬ 
tions  for  varying  number  densities  if  desired. 
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6 .  RADIATIVE  TRANSFER 
6.1  Technique  Selection  Criteria 

Radiative  transfer  theory  provides  a  mathematical 
description  of  the  interaction  between  incident  solar  radia¬ 
tion  and  the  relevant  optically  active  constituents  of  the 
atmosphere.  For  our  purposes,  these  species  include  both 
ambient  gases  and  aerosols.  The  interaction  mechanisms 
treated  by  radiative  transfer  theory  have  been  alluded  to 
in  the  predeeding  discussions  of  atmospheric  optical  proper¬ 
ties  and  include  molecular  (Rayleigh)  scattering  (§3.1), 
gaseous  absorption  (§3.2),  and  aerosol  absorption  and 
scattering  (§5.2).  For  the  purpose  of  this  study,  it  is 
minimally  required  that  the  radiative  transfer  treatment 
adopted  handle:  (a)  multiple  scattering,  (b)  the  inherent 
anisotropic (i . e .  highly  directional)  scattering  characteristic 
of  aerosols,  (c)  reflection  of  radiation  at  the  atmosphere- 
ocean  interface,  and  (d)  the  azimuthal  dependence  given  by 
the  sun/sensor/geometry. 

A  hierarchy  of  potential  radiative  transfer  treatments 
exists  (Hansen  and  Travis,  1974;  Lenoble,  1977)  ranging  from 
the  simple,  single  scattering  analysis  demonstrated  in  Fett 
and  Isaacs  (1979)  to  a  variety  of  highly  numerical  computa¬ 
tional  algorithms  (e.g.  Braslau  and  Dave,  1947  a,b).  In 
selecting  an  appropriate  approach  in  this  investigation, 
criteria  have  included:  (a)  ability  to  satisfactorily 
treat  mechanisms  cited  above,  (b)  maximization  of  the  number 
of  sensitivity  analyses  which  could  be  undertaken,  (c) 
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optimization  of  available  computer  resources,  and  (d) 
consistency  of  approach  with  accuracies  in  other  components 
of  the  investigation  such  as  in  specifying  the  physical 
aerosol  models.  For  these  reasons  an  approximate,  analytical 
treatment  was  adopted. 

6.2  Analytical  Approaches 

From  a  practical  perspective,  the  essential  difference 
between  numerical  and  analytical  approaches  is  that  the 
latter  require  very  little  computational  effort  and  hence  a 
significant  amount  of  computer  time  may  be  saved  if  they  are 
implemented.  Unfortunately,  exact  analytical  treatments  are 
available  only  for  a  few  cases  which  are  not  immediately 
applicable  to  geophysical  remote  sensing  problems.  They  are 
of  interest  mathematically,  however,  to  compare  with 
corresponding  cases  of  more  computationally  complex  numerical 
techniques.  These  cases  include,  for  the  most  part,  approaches 
based  on  Chandrasekhar's  H  functions  for  semi-infinite 
atmospheres  and  X  and  Y  functions  for  finite  atmospheres 
(Chandrasekhar,  1950) .  For  the  former  set  of  problems 
solutions  are  available  for  isotropic  (Stibbs  and  Weir,  1959; 
Abhyankar  and  Fymat,  1971),  Rayleigh  (Chandrasekhar,  1950; 
Lenoble,  1970).  and  various  anistropic  phase  funstions 
(Chandrasekhar,  1950;  Sobolev,  1956;  Kolesov,  1972).  For 
finite  atmospheres,  solutions  are  available  only  for  isotropic 
(Carlstedt,  1966)  and  Rayleigh  scattering  (Sekera  and  Kahle, 
1966)  . 
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For  remote  sensing  problems  relevant  to  radiative 
transfer  in  the  atmosphere-ocean  system,  methods  are  re¬ 
quired  which  treat  a  finite  atmosphere  (i.e.,  the  total 
optical  depth  is  not  infinite)  with  scattering  properties 
which  are  anisotropic  in  order  to  simulate  aerosol 
scattering.  Although  exact  analytical  treatments  are  not 
available  in  such  cases,  a  variety  of  approximate  analytical 
methods  with  quantifiable  accuracy  may  be  employed.  The 
utility  of  approximate  analytical  treatments  lies  in  their 
extreme  computational  efficiency  while  retaining  many  of  the 
salient  physical  mechanisms  of  radiative  transfer  involved. 

Approximate  analytical  methods  may  be  conveniently 
classified  into  one  of  two  groups: 

.  approaches  based  on  taking  the  first  few  tractable 
orders  of  more  extensive  numerical  treatments 

.  approaches  formulated  specifically  as  approximate 
treatments 

Examples  falling  into  the  first  category  include  first 
(Deirmendjian,  1969)  and  second  (Hovenier,  1971)  order  of 
scattering  treatments  explicitly  formulated  from  successive 
order  of  scattering  approaches  (Irvine,  1965;  Nagel  et  al., 
1978),  analytic  two  and  four  steam  (Liou,  1974)  approxima¬ 
tions  based  on  discrete  ordinate  methods  (Liou,  1973) ,  and 
the  two-step  function  approach  (Burke  and  Sze,  1977)  deriv¬ 
ed  from  more  general  variational  methods  (Sze,  1976). 

Treatments  explicitly  formulated  as  approximate 
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approaches  include  various  similarity  relations  (van  de  Hulst 
and  Grossman,  1968;  Hansen,  1969),  the  Eddington  approxima¬ 
tion  (Shettle  and  Weinman,  1970),  and  various  general  two- 
stream  methods  (Chu  and  Churchill,  1955;  Coakley  and 
Chylek ,  1975). 

The  degree  of  accuracy  intrinsic  to  a  specific  approx¬ 
imate  analytical  treatment  varies  not  only  with  the  treat¬ 
ment  itself,  but  with  the  relevant  optical  propagation  par¬ 
ameters  involved.  Expected  errors  may  be  quantified  by  ex¬ 
amining  certain  standard  cases  and  comparing  accuracies 
either  with  available  exact  solutions  or  with  numerical" 
solutions  of  specified  precision.  For  example.  Table  9 
lists  percent  error  figures  for  a  comparison  between 
emergent,  I(0,-1.0,y),  and  transmitted,  I(t*,  -1.0,y), 
intensities  derived  from  single  scattering  vs.  multiple 
scattering  (Coulson  et  al. ,  1960)  treatments  in  a  Rayleigh 
atmosphere  with  the  sun  at  zenith  (yo=1.0)  and  zero  surface 
albedo  (AL=0.0).*  Results  are  presented  for  various  observer 
zenith  angles  (y)  and  total  atmospheric  optical  depths  ( x* ) . 
Note  that  even  for  this  fundamental  approximate  analytical 
approach,  errors  are  highly  dependent  on  atmospheric  optical 
depth,  emphasizing  the  need  to  quantify  the  behavior  of 
such  treatments  a  priori. 


These  parameters  are  defined  in  §6.4. 
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Table  9 
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(T  =  0 

/M) 
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/  * 
(T  =  T 
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T* 

★ 
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0.10 

0.25 

0.50 

0.10 

0.25 

0.50 

M 

1.00 

13.2 

24.6 

37.1 

12.3 

24.7 

37.9 

0.72 

13.8 

25.6 

39.3 

13.4 

26.1 

40.3 

0.52 

14.3 

26.9 

40.8 

13.9 

27.2 

42.3 

0.28 

14.7 

27.8 

41.9 

14.6 

28.5 

44.6 

0.10 

14.8 

27.6 

40.2 

14.9 

29.3 

46.1 

Percent  error:  single  scattering  vs.  multiple  scattering 
with  Rayleigh  phase  function  (y„  =  1.0,  A  =  0.0). 
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6 . 3  Adopted  Approach 

Of  the  approximate  techniques  described  above,  simple 
finite  stream  approximations  are  particularly  convenient  to 
use.  A  recent  review  (Meador  and  Weaver,  1980)  of  extant 
two-stream  approximations  has  compared  the  many  approaches 
available  within  the  formalism  and  defined  domains  of 
applicability  and  accuracy  vs.  total  optical  depth.  Two- 
stream  approximations  are  attractive,  since  accuracies  are 
easily  quantified,  and  the  mathematics  reduces  to  a  set  of 
analytical  (although  algebraically  complex)  equations.  For 
the  latter  reason,  computational  efficiences  are  extremely 
high  (i.e.,  machine-usage  time  is  minimal).  Furthermore, 
with  appropriate  care  in  the  formulation,  scattering  aniso- 
topy  and  surface  reflection  may  be  treated  by  employing 
physically  realistic  boundary  conditions. 

For  these  reasons,  the  specific  treatment  adopted  in 
this  work  is  based  on  a  recent  extension  (Kaufman,  1979)  to 
an  earlier  two-stream  formulation  (Coakley  and  Chylek,  1975) , 
which  allows  the  specific  azimuth/zenith  angle  dependence  of 
emergent  and  transmitted  intensities  to  be  simulated.  By 
employing  calculated  backscatter  fractions  (Wiscombe  and 
Grams,  1976)  for  general  aerosol-phase  functions  in  the 
analysis,  the  general  scattering  anisotropy  of  the  real 
atmosphere-ocean  system  may  be  calculated.  This  particular 
two-stream  approach  (Coakley  and  Chylek,  1975)  is  formulated 
to  be  well-behaved  in  the  thin  atmosphere  limit,  a  trait 
particularly  valuable  for  treating  problems  in  the  visible 
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and  near  infrared  spectral  regions. 


6 . 4  Theory 

In  order  to  simulate  the  emergent  radiance  I , 
measured  by  a  meteorological  satellite  viewing  the  earth- 
atmosphere  system,  it  is  necessary  to  solve  the  radiative 
transfer  equation  for  each  wavelength  (Goody,  1964) :+ 

cl  I 

y  (T,y,<|>)  =  I  (t,u,<M-u>0  Js(t,u,4>)  (33) 

where : 

Js  ("r  ,y  ,<J>)  = 

-t/u  i  r2Tr+fi 

- -  Fe  /Mo  +  J  J  P(e)I(T,y;*)dy'd*'  (34) 

0  -1 

In  the  above: 

_  n 

I  (r  ,y,4>)  =  wavelength  dependent  radiance  (mw  cm 

ym  ^  sr  ^) 

J(t  ,y,4>)  =  wavelength  dependent  source  function 

CO 

*=£'  6^(z)dz  =  wavelength  dependent  optical  depth 

it  co 

t=  /  $^(z)dZ  =  total  optical  depth  from  surface  to  space 
o 

+  ...  » 

The  derivation  is  treated  in  greater  detail  in  the  Appendix. 


8^  (z)  =  6®  (X)  +  8g(X)  +  0e(X)  =  height  dependent, 

total  extinction  coefficient  per  unit  length 
(km  including  Rayleigh  scattering,  gas 
absorption,  and  aerosol  extinction  contribu¬ 
tions  . 

U  =  cos0  (0  zenith  angle  of  emergent  radiance) 
y0  =  cos0o  (0O  zenith  angle  of  sun) 

<j>  =  azimuth  angle  of  emergent  radiance 
w0  (X)  =  wavelength  dependent  single  scattering  albedo 
P  [0]  =  wavelength  dependent  phase  function  for 
redirection  of  scattered  radiance^ 

0  =  scattering  angle 

=  yy'  +  (1-y2)  **  (l-y'^^cos  (<f>— 4> ' ) 
i.e.  incident  beam  y,<j> 

scattered  beam  y,<J>‘ 

0 o  =  scattering  angle  between  direction  of  emergent 
radiance  and  direction  of  solar  incidence 
(yo •  ^ 

1The  phase  function  P(0)  is  the  effective  phase  function  due 
to  both  molecular  (Rayleigh)  scattering  and  aerosol  scatter¬ 
ing  properly  weighted: 

P  (0)  -  b|(z)Pr(0)  +  3s(z)Pa(0) 

8|(z)  +  8s (z) 

where 

6  s  g  s 

R,p  =  scattering  coefficients  due  to  Rayleigh  and 
aerosol  scattering,  respectively 

PR(0).PA(0)=  phase  functions  for  Rayleigh  and  aerosol  compon¬ 
ents 


F  =  wavelength  dependent  extraterrestrial  solar 
flux 

(These  values  shown  in  Fig.  21  from  Thekaekara 
et  al,  1969) 

Equations  (33)  and  (34)  may  be  transformed  into  coupled 
differential  equations  for  upward,  F+(t),  and  downward, 

F  (t) ,  flux  profiles: 


AF  =  f 


where  F  = 


rF-(T)  i 

lF  (t)  j 


f  = 


[fl]  , 


(35) 


the  operator  A  is  given  by: 


A  = 


(36) 


where  1  is  the  identity  operator: 


Y1  =  2  [1  “  U)o(1  "  6,)  1  ' 

Y  2  =  2u)q  8  1  , 

fl  =  -TTFw0f3(y0)e-T/yo  , 

f2  =  ttFu)0[1  -  8(p0)  ]e"T/y° 


(37) 


and 


8'  =  backscatter  fraction  for  isotropically 

incident  radiation  (Wiscombe  and  Grams,  1976) 
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Wavelength  (pm) 

Figure  21 

Solar  Flux  (Thekaekara  et  al.,  1969) 
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B(m)  backscatter  fractions  for  raonodirectionally 

3(y0)  incident  radiation  (Wiscombe  and  Grams,  1976) 

The  backscatter  fractions  defined  above  result  from 
angular  integration  over  the  phase  function  in  Equation  (34) 
and  are  dependent  on  the  asymmetry  parameter  (g)  for  the 
specific  wavelength  under  consideration.  Figure  22 
illustrates  the  dependence  of  both  6'  and  3(y)  on  asymmetry 
parameter.  Generally  3'  varies  from  0.5  for  isotropic  or 
Rayleigh  scattering  (g=0)  to  0.0  for  total  forward 
scattering  (g=1.0)  [i.e.  there  is  no  backscatter.)  The 
range  of  6(y)  values  is  similar,  however,  the  rate  of 
charge  in  the  domain  0 . 0  <_  g  £  1 . 0  is  dependent  on  the  speci¬ 
fic  value  of  y  (=  cos0)  chosen. 

Equation  (35)  is  subject  to  the  boundary  conditions: 

F~(t  =  0)  =  0 

F+  ( x  =  x*)  =  AL(TTFy0  e_T*/y°  (38) 

+  F~ (t  =  t*)  ] . 

where  AL  is  the  wavelength-dependent  surface  albedo.  For  the 
purpose  of  this  work,  the  surface  has  been  assumed  to  be  a 
Lambert  reflecter  with  an  albedo  spectrum  taken  from  Ramsey 
(1968) [as  given  in  Curran  (1972)].  Ocean  surface  albedo  as 
a  function  of  wavelength  is  illustrated  in  Figure  23. 
Solutions  to  (35)  have  the  simple  forms: 

F  =  A  exp(kx)  +  B  exp  (-kx)  +  C  exp  (-x/y0)  (39) 

where  k  represents  the  eigenvalues  of  (35) . 
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Asymmetry  Parameter,  g 


Figure  22 

Backscatter  Fractions  as  a  Function 
of  Asymmetry  Parameter 
(Wiscombe  and  Grams,  1976) 
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(Ramsey,  1968) 
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The  zenith  angle  dependent  radiances  I  (T,y,<j>)  will 
satisfy : 


LI  =  -  —  CF  -  f-  e_T/yo  p 

~  it  ~  4 


(40) 


where  I  = 


F  = 


(t  ,u,4>) 
LI  (x ,y ,<j>) 


[F-  (T)1 
lF  (t)  j 


p  =  rF(u/— yQ/d)i 
LP(y,y0,<j>)  J 

and  the  operator  L  is  given  by: 


and 


*  =  ^  dT  -  il 


r  =  r [1-6 (u) ]  8<y) , 

~  l3(y)  [l-(3(y)]J  • 


(41) 


(42) 


Equation  (40)  is  subject  to  boundary  conditions  derived 
from  (38)  namely: 


I  (x  =  0  ,y ,  (f>)  =  0.0 


(43) 


and 


I+(x  =  x*,y,<i>)  =  F+(t  =  t*)/tt 


(44) 


Solutions  to  (40)  have  the  simple  form: 

!  =  d  exp  (  +  x/y)  +  V  exp(-kT) 

u  exp^-x/yj  +  (1  _  yk)  +  (1  +  yk) 


W  exp (-T  y0) 


(1 


P  V  t  MQ) 

+  m/m0) 


(45) 
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where  D,  U,  V,  W  are  analytical,  algebraic  functions  of  the 
input  variables. 

A  variety  of  comparisons  with  exact  and  numerical 
results  have  been  made  to  quantify  the  accuracy  of  this 
simple  analytical  approach.  Table  10  summarizes  a  comparison 
between  the  two-steam  model  and  exact  results  for  isotropic 
scattering  obtainable  using  Chandrasekhar's  X  and  Y  functions 
tabulated  in  Carlstedt  and  Mullikin  (1966) .  It  is  particu¬ 
larly  notable  that  errors  are  less  than  5%  over  much  of  this 
domain,  and  particularly  for  emergent  intensities  with  zenith 
angles  approaching  unity  (y+1.0).  This  geometry  simulates  a 
nadir-pointing  satellite. 

Analogous  results  for  Rayleigh  scattering  are  provided 
in  Table  11.  In  this  case,  the  solar  zenith  angle  is  fixed 
at  57°  and  the  observer  zenith  at  13.5°  to  simulate  a 
satellite  field  of  view  for  a  polar-orbiting  sensor.  Percent 
errors  are  given,  comparing  two-stream  results  to  the  exact 
calculation  of  Rayleigh  scattering  by  Coulson  et  al.  (1960) 
for  a  variety  of  optical  depths,  surface  albedos  of  0.0  and 
0.25,  and  azimuth  distances  of  30°,  90°,  and  150°.  A 
comparison  of  these  results  to  the  single  scattering  approxi¬ 
mation  results  in  Table  9,  indicates  that  a  much  higher  degree 
of  accuracy  is  achievable  using  the  two-stream  approach. 

The  most  stringent  test  of  an  approximate  multiple¬ 
scattering  radiative  transfer  model  is  highly  anistropic 
aerosol  scattering.  Results  of  comparison  with  numerical 
successive  order  of  scattering  approaches  (Nagel  et  al.,1978) 
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Table  10 
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Table  11 
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essentially  verify  (Kaufman,  1979)  accuracy  on  the  order  of 
6-18%  for  the  near  zenith  observation  angles  characteristic 
of  many  polar-orbiting  satellites.  Accuracy  for  these  cases 
degrades  to  unacceptable  levels  for  solar  or  observer  zenith 
angles  near  the  horizon.  However,  for  non-terminator  solar 
illumination  conditions,  the  two-stream  approach  should 
provide  a  practically  useful  remote  sensing  analysis  tool. 


6.5  Ensemble  Atmosphere  Optical  Properties 

The  total  optical  depth,  x*,  single  scatter  albedo,  w0 , 
and  scattering  asymmetry  parameter,  g,  (which  determines  the 
backscatter  fractions,  see  Figure  22)  used  in  solving 
equations  (35)  to  (45)  are  derived  from  the  wavelength  and 
meteorologically  dependent  atmospheric  optical  properties 
discussed  previously  (see  Figure  2) .  These  ensemble 
atmospheric  optical  properties  are  given  by: 


x*(A)  =  x|(A)  +  x  (A)  +  x®(A)  +  r®(A) 
w0(A)  =  x|(A)  +  x®(A)  +x®(A)  /x*(A) 


g  (A)  = 


i  *  ( A  )  Cjx  (A)  +  Tf*(A)g2(A) 
T®(A)  +  (A)  +  (A) 


where  g^  and  g2  are  asymmetry  parameters  for  the  accumulation 
and  coarse  modes,  respectively.  The  optical  depths  are  those  for 

g 

Rayleigh  scattering  [iR;  eqn .( 4 )], gaseous  absorption[x  ;eqn. (6) ] , 

0  g 

and  aerosol  extinction  (x^  and  scattering  (x^  2)  for  the  accum¬ 
ulation  and  coarse  modes  respectively  (eqn.  (32)].  Eqn.  (48) 


represents  a  weighting  of  the  respective  Rayleigh  and 
aerosol  scattering  phase  function  with  g (Rayleigh)  =  0.0. 


6.6  DMSP  Bandpass  Weighted  Radiances 

The  solution  to  equation  (45)  for  each  set  of  meteoro¬ 
logical  variables  investigated  provides  a  wavelength- 
dependent  intensity  spectrum,  I  (X) .  The  intensity 
measured  by  the  DMSP  sensor  at  the  satellite,  I  (DMSP) , 
however,  is  given  by  weighting  these  monochromatic  intensi¬ 
ties  by  the  DMSP  spectral  bandpass  function  given  in  Figure  1 
according  to: 


I  (DMSP) 


4>  ( A )  I  ( A )  dA 


AA 


(49) 


where  4>  is  the  sensor  response  function  over  wavelength 
interval  AA  (0. 4-1.1  um  for  DMSP  VHR,  LF) .  (See  Figure  1.) 
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7.  DISCUSSION  OF  RESULTS 


7.1  Scope 

In  support  of  the  objectives  of  this  study,  the 
sensitivity  of  DMSP  incident  radiances  in  the  0. 4-1.1  ym 
spectral  region  to  variations  in  meteorological  variables 
(see  Table  1)  was  investigated  by  performing  approximately 
3600  individual  calculations  employing  the  technical 
approach  summarized  in  Figure  2  and  described  in  Sections  3 
to  6.  These  results  include  both  wavelength-dependent, 
simulated  radiance  spectra,  I (A),  and  DMSP-weighted 
radiance,  I  (DMSP).  A  discussion  of  these  results  follows. 

7.2  Maritime  Model  Geometry 

As  discussed  in  §6.4  above,  radiance  computations  are 
dependent  on  the  relative  sun/sensor  geometry  through 
various  backscatter  fractions  and  the  composite  angular 
scattering  (phase)  function  of  the  atmosphere.  For  the 
purpose  of  this  study  these  geometric  factors  were 
deemuhasized  in  favor  of  the  meteorological  dependences 
of  the  results.  Figures  24  and  25,  however,  illustrate  the 
dependence  of  radiances  on  solar  zenith  angle  cosine 
(=  cos  ’)  )  and  relative  azimuth  angle,  respectively.  The 
model  atmosphere  is  the  maritime  model  with  RH  =  80%  and  a 
scale  height  of  1.2  km.  A  nadir  viewing  sensor  (y=1.0)  is 
assumed.  Results  in  Figure  24  indicate  a  general  increase 
in  radiances  when  the  sun  is  closer  to  the  zenith  position 
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(pQ=1.0)+.  Since  this  effect  is  dependent  on  optical  depth, 
it  is  enhanced  toward  shorter  wavelengths  where  molecular 
scattering  is  important.  Since  most  of  the  available  solar 
radiation  is  at  shorter  wavelengths  (see  Figure  21) ,  DMSP 
radiances  are  increased  by  this  effect. 

Since  the  general  spectral  features  of  the  radiances 
is  Figure  24  are  not  markedly  different  for  the  two  cases, 
subsequent  calculations  assume  pq=0.5. 

Radiances  for  a  nadir  viewing  instrument  (p=1.0)  are 
independent  of  azimuth  angle  difference  [  see 

definitions  related  to  equations  (33)  and  (34)].  For  off- 
nadir  (p.^1.0)  positions  within  the  sensor  scan,  radiances 
may  vary  with  azimuth  angle.  In  Figure  25,  a  sensor  zenith 
angle  of  about  37°  (p=0.8)  is  assumed  and  azimuth  angles  of 
0,45,90,135,  and  180°  are  evaluated  for  a  solar  zenith  angle 
cosine  of  0.8.  The  five  curves  are  generally  coincident 
except  for  a  small  region  in  the  vicinity  of  0.4  pm.  These 
calculations  indicate  that  DMSP  radiances  are  approximately 
independent  of  azimuth  angle  difference. 

7.3  Maritime  Model  Surface  Albedo 

The  Lambert  surface  albedo  assumption  adopted  in  the 
radiative  transfer  theory  boundary  conditions  [  equation 
(38)  ]  ignores  the  possibilities  of  surface  albedo  enhancement 

-t-Note:  The  solar  zenith  angles  in  Figure  24  are  0  =cos  ^p 

or  60°  and  25.8°  for  po  =  0.5  and  0.9,  respectively. 
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due  to  sun  glitter  on  a  wind-ruffled  ocean  (Guinn  et  al, 
1979)  or  decreases  due  to  subsurface  absorption  by  biologi¬ 
cal  materials  such  as  phytoploukton  (Curran,  1972)  .  These 
effects  were  parametrically  investigated  by  varying  the 
assumed  ocean  surface  albedo  (see  Figure  23)  with  multipli¬ 
cative  factors  of  0.0,  0.5,  1.0,  1.5,  and  2.0,  where  1.0 
represent  the  assumed  surface  albedo  and  2.0  an  increase  by 
a  factor  of  2.0.  These  results  are  illustrated  in  Figure 
25(a).  It  is  noteworthy  that  any  variations  in  surface 
albedo  affect  only  the  visible  portion  (X<0.7ym)  of  DMSP 
incident  radiances. 

Increased  surface  albedos  increase  radiances  from  0.4- 
0.7  urn  [Figure  25(b)]  although  the  effect  is  not  linearly 
dependent  on  the  surface  albedo.  Bandpass-weighted  calcu¬ 
lations  corresponding  to  those  in  Figures  25(a)  and  (b) 
indicate  a  20%  enhancement  in  radiance  for  a  change  from 
zero  albedo  to  2.0  times  the  assumed  value. 

7.4  Maritime  Model  -  Scale  Height  and  Relative  Humidity 

Dependence 

A  fundamental  objective  of  this  study  was  to  present 
an  internally  consistent  picture  of  simulated  radiance 
dependence  on  relative  humidity  including  the  effects  of 
both  water  vapor  absorption  and  aerosol  growth.  Figures 
26(a)  -  (d)  present  the  results  of  calculations  illustrating 
the  dependence  of  the  maritime  aerosol  model  (N=4000  cm  ^) 
on  relative  humidity  for  assumed  scale  heights  of  0.5, 1.0, 
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1.5,  and  2.0  km,  respectively+ .  Each  figure  consists  of  five 
curves  corresponding  to  50,70,80,90,  and  95%  relative  humidity. 

For  the  fixed  number  density  above,  surface  visual  ranges  in 

the  modeled  situations  vary  with  relative  humidity  according 

to  Table  7 . 

An  apparent  feature  of  these  results  is  the  marked 
increase  of  simulated  radiances  in  the  near  infrared  window 
regions  (see  §3.2)  with  increases  in  relative  humidity.  In 
contrast,  this  response  is  much  less  evident  at  visible  wavelen¬ 
gths.  For  example,  assuming  a  scale  height  of  0.5  [Figure 
26 (a) ] ,  the  ratio  of  simulated  intensities  for  relative 
humidities  of  95  (visual  range,  8.7  km)  and  50%  (visual 
range,  39.6  km)  is  near  unity  at  0.55  pm,  while  it  is  about 
2.8  near  1.0  pm.  This  corresponds  to  an  approximate  28% 
increase  in  bandpass-weighted  DMSP  radiance  for  the  transi¬ 
tion  from  50  to  95%  relative  humidity.  These  results  are  in 
qualitative  agreement  with  those  discussed  in  Fett  and 
Isaacs  (1979).  As  scale  height  increases,  there  is  a 
corresponding  increase  in  radiances  throughout  the  DMSP 
spectral  interval.  For  a  fixed  relative  humidity  of  80%, 
for  example,  a  scale  height  increase  from  0.5  to  2.0  km 
corresponds  to  an  increase  in  DMSP  bandpass-weighted 
radiances  by  a  factor  of  about  2.0. 


+Note:  As  discussed  in  §5.4,  variations  in  scale  height  may 

also  be  interpreted  as  increase  in  number  density 
(with  corresponding  charges  visual  range)  for  a  fixed 
scale  height. 
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In  light  of  the  strong  dependence  of  visual  range  on 
relative  humidity-dependent  aerosol  extinction  as  given  by 
the  Koschmieder  relation  [Equation  (30)]  applied  in  Table  7, 
the  apparent  insensitivity  of  simulated  emergent  visible 
radiances  to  changes  in  relative  hunidity  may  be  difficult 
to  comprehend-  The  situation  may  be  partially  reconciled 
by  examination  of  the  back  scattering  geometry  considered  in 
the  radiative  transfer  calculation  and  the  angular  scatterinq 
properties  of  the  aerosols  treated  in  comparison  to  Rayleigh 
scattering  by  gases.  While  visual  range  is  inversely  propor¬ 
tional  to  the  sum  of  Rayleigh  scattering  and  aerosol  extinc¬ 
tion  coefficients  visible  radiances  are  qualitatively  depend¬ 
ent  on  the  product  of  this  sum  and  the  appropriate  angular 
scattering  function.  As  described  in  Section  5,  increasing 
relative  humidity  implies  not  only  increased  aerosol  extinc¬ 
tion  [see  Figure  15(a)],  but  also  larger  asymmetry  parameters 
[see  Figure  15(b) ].  Larger  asymmetry  parameters  mean  less 
back  scattering  to  the  satellite  sensor.  At  visible  wave¬ 
lengths,  the  effective  angular  scattering  function  (see 
footnote  p.  6-8)  or  asymmetry  factor  [Equation  (48)]  must  be 
weighted  to  account  for  both  aerosol  scattering  and  Rayleigh 
scattering.  Rayleigh  scattering  back  scatters  much  more  than 
aerosol  scattering  at  the  scattering  angles  of  120°  used  here. 
Therefore,  although  increases  in  relative  humidity  increase 
aerosol  extinction,  they  decrease  the  effective  magnitude  of 
backscattering  due  to  aerosols  and  gases.  These  effects 
partially  compensate  for  one  another  in  the  visible  region. 
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At  near  infrared  wavelengths,  however,  the  influence  of 


Rayleigh  scattering  is  small.  Thus,  the  backscatter  is 
primarily  determined  by  the  aerosol  angular  scattering 
function  and  increases  in  aerosol  extinction  are  directly 
translated  into  increased  radiances  without  compensatory 
decreases  in  the  effective  amount  of  backscatter 
7.5  Hybrid  Model  -  Relative  Humidity,  Wind  Speed,  and  Size 

Range  Dependence 

As  discussed  in  §4.3  the  hybrid  tropospheric-Munn  Katz 
model  provides  the  capabilities  to  examine  dependence  of 
emergent  radiances  on  size  range  and  wind  speed  in  addition 
to  scale  height  and  relative  humidity. 

7.5.1  Hybrid  Model  -  Accumulation  Mode  only 

One  question  which  arose  early  in  this  study  concerns 
the  relative  role  of  fine  (accumulation)  mode  and  coarse 
mode  aerosols  in  determining  wavelength-dependent  emergent 
radiances.  In  order  to  address  this  point,  simulations  were 
performed  using  only  the  accumulation  mode  of  the  adopted 
hybrid  model  (N=5000cm  ^).  These  results  are  presented  in 
Figures  27 (u) -(d)  as  functions  of  relative  humidity  (50,70, 
80,90,95?,)  for  scale  heiyhts  of  0.5,  1.0,  1.5,  and  2.0  km. 
Visual  ranges  corresponding  to  these  relacive  humidities  are 
given  in  Table  7.  Notable  in  these  figures  is  the  suppressed 
dependence  of  accumulation  mode  aerosols  on  relative 
humidity  when  compared  to  the  bimodal  maritime  results 
presented  in  57.4.  This  is  particularly  true  in  the  near 
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infrared  window  regions.  A  comparison  between  corresponding 
cases  illustrated  in  Figures  26  and  27  suggests  that  coarse 
mode  aerosols  play  a  predominant  role  in  influencing  the 
dependence  of  simulated  emergent  radiances  on  relative 
humidity.  For  moderate  relative  humidities  (i. e. <70-80%) , 
visible  radiances  are  much  more  dependent  on  accumulation 

mode  aerosols  while  near  infrared  wavelengths  are  influenced 
primarily  by  the  coarse  mode. 

7.5.2  Hybrid  Model  -  Relative  Humidity 

Figures  28  (a) -(d)  illustrate  the  effect  of  including 
the  coarse  mode  of  the  hybrid  aerosol  model.  The  accumula¬ 
tion  mode  is  as  described  above  with  number  densities  in  the 
coarse  mode  determined  by  relative  humidity  and  wind  speed  as 
given  in  Table  6.  The  visual  ranges  corresponding  to  these 
meteorological  variable  for  the  combined  model  are  given  in 
Table  7.  For  these  cases  scale  heights  are  fixed  at  1.2  km 
(HI)  for  the  accumulation  mode  and  0.8  km  (H2)  for  the 
coarse  mode.  Suitable  comparisons  (H=1.0)  are  to  Figure 
26 (b)  for  the  maritime  model  or  Figure  27 (b)  for  the 
accumulation  mode  only  calculation. 

Dependence  is  shown  as  a  function  of  relative  humidity 
(50,70,80,90,95%)  for  fixed  windspeeds  of  0,5,7,  and  10  ms-1, 
respectively.  The  influence  of  coarse  mode  aerosols  on  near 
1R  radiances  for  this  model  is  approximately  intermediate 
between  those  for  either  the  maritime  model  (which  generally 
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Log [Radiance] 


Figure  28 (c) 

Dependence  on  Relative  Humidity 
(Hybrid  Model,  Hi  =  1.2  km, 

H2  =  0.8  km,  WS  =  7  ms-1) 


Figure  28(d) 

Dependence  on  Relative  Humidity 
(Hybrid  Model,  Hi  =  1.2  tan, 

H2  =  0.8  kn,  VIS  »  10  no"1) 
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has  a  higher  coarse  mode  number  density)  or  the  tropospheric 
model  (which  has  no  coarse  mode) .  Again,  there  is  very 
little  dependence  of  visible  radiances  on  relative  humidity. 

7.5.3  Hybrid  Model  Wind  Speed 

For  fixed  relative  humidities,  increases  in  wind  speed 
enhance  simulated  radiances  particularly  in  the  near  IR 

region.  Figure  29 (a) -(c)  illustrate  the  dependence  of 
radiance  on  windspeed  for  relative  humidities  of  50,  80,  and 
95%,  respectively.  Again,  the  effect  of  windspeed  is  not 
apparent  at  visible  wavelengths.  Figures  28  and  29 
collectively  suggest  that  both  windspeed  and  relative 
humidity  are  crucial  in  determining  near  IR  radiances 
incident  on  the  DMSP  censor. 

7.6  Visual  Range 

The  discussion  in  §5.3  provides  a  prescription  to 
adjust  number  densities  in  the  adopted  models  to  correspond 
to  specific  surface  visual  ranges  at  fixed  relative  humidity 
(and  if  applicable  wind  speed) .  These  number  densities  are 
summarized  in  Table  8.  Calculations  were  performed  using 
these  values  at  visual  ranges  of  5,10,23,  and  50  km  for  both 
the  maritime  and  hybrid  models.  These  results  are  illustra¬ 
ted  in  figures  30  and  31,  respectively.  Results  are  quali¬ 
tatively  similar,  indicating  an  increase  in  radiances  at  all 
wavelengths  between  0.4  and  1.1  ym  with  a  decrease  in  visual 
range.  Note  that  visible  wavelengths  are  dramatically 


Log [Radiance]  Log [Radiance]  Log [Radiance] 


Wavelength  (pm) 


(a)  RH  =  50% 

i  ► 


I  V  s  V  s  U  S 


Wavelength  (yin) 


(b)  RH  =  80% 


‘i  .6  8  ' 

Wavelength  (pm) 


(c)  RH  -  95% 

Figure  29 

Dependence  on  Wind  Speed 
(Hybrid  Model,  Hi- 1.2  km,  II^O.o  km) 
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Figure  30 

Dependence  on  Visual  Range 
(Maritime  Model/  II  =  1.2  km,  RH  =  80%) 


affected  since  the  number  density  adjustments  cited  above 
are  applied  to  both  the  accumulation  and  coarse  modes  in  the 
maritime  model  and  exclusively  to  the  fine  mode  in  the  hybrid 
model . 

7.7  DMSP  Bandpass-Weighted  Radiances 

In  the  preceedinq  discussion  of  results  reference  was 
made  to  DMSP  bandpass-weighted  radiances.  For  each  of  the 
cases  described  above,  simulated  DMSP  radiances  were 
evaluated  according  to  the  procedure  described  in  §6.6. 
Plotting  these  values  on  a  common  scale,  the  sensitivity 
of  simulated  DMSP  radiances  to  variations  in  the  meteorolo¬ 
gical  variables  discussed  above  can  be  demonstrated.  These 
results,  although  preliminary,  also  provide  potentially  useful 
tools  for  meteorological  analyses  which  may  be  developed 
into  user  oriented  operational  nomographs. 

7.7.1  Maritime  Model 

Figure  32,  illustrates  the  dependence  of  DMSP  bandpass- 
weighted  radiances  on  both  relative  humidity  (RH)  and 
scale  height  (H)  for  the  maritime  model.  As  previously 
discussed  (§5.4),  scale  height  may  act  as  a  surrogate  for 
changes  in  number  density  with  a  fixed  scale  height.  Solid 
curves  represent  the  variation  of  radiance  with  scale  height 
(upper  ordinate)  for  the  labeled,  fixed  relative  humidity. 
Dashed  curves  present  the  converse  case,  that  is,  variation 
with  relative  humidity  (lower  ordinate)  for  the  labeled, 
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Bandpass-Weighted  Radiances  (Maritime  Model) 


fixed  scale  height.  As  indicated  in  the  figure,  a  particular 
DMSP  radiance  is  a  multivalued  function  of  relative  humidity 
and  scale  height.  For  example,  assuming  I  (DMSP)  =  1.0  mWcm 
vim'^sr-^ ,  this  radiance  corresponds  to  (RH,H)  pairs  including 
(95,1.0),  (90,1.3),  (80,1.8),  etc.  Therefore,  it  is  unlikely 

that  both  values  can  be  uniquely  determined  from  DMSP  data 

alone.  However,  supporting  meteorological  data  or  other 

remotely  sensed  data  may  provide  information  on  either 

variable.  In  this  case,  DMSP  radiances  can  be  used  to 

—  -2 

estimate  the  other.  For  example,  if  I  (DMSP)  =  1.0  mWcm 

ym  ^sr  ^  and  RH=80%,  the  scale  height  is  about  1.8  km.+ 

Furthermore,  if  it  is  estimated  that  the  actual  scale  height 

is  1.2  km,  then  there  must  be  about  1.5  (=1.8/1. 2)  times 

more  aerosol  particles  present  than  in  the  aerosol  model 
-3  -3 

(i.e.,  4000  cm  )  or  6000  cm  .  For  80%  relative  humidity. 
Table  8  suggests  a  surface  visual  range  of  slightly  over 
10  km. 


7.7.2  Hybrid  Model 

The  introduction  of  wind  speed  dependence  complicates 

the  picture  somewhat  for  the  hybrid  model,  although  the  basic 

idea  remains  the  same.  Figure  33,  illustrates  these  results 

with  labeled  curves  corresponding  to  fine  mode  scale  height 

(Hi),  the  coarse  mode  scale  height  is  fixed  at  0.8  km, 

+Note:  Due  to  the  approximate  nature  of  the  model  used  in 
these  calculations,  such  applications  will  have 
computational  errors  associated  with  them. 
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dpass- 


relative  humidity  (RH)  ,  and  wind  speed  (VJS)  .  For  fixed 

values  of  two  of  these  variables,  the  appropriate  estimate  of 

the  third  corresponding  to  a  specific  DMSP  radiance  can  be 

read  from  the  relevant  ordinate  scale.  For  example,  if  T 

-2  -1  -1 

(DMSP)  =  1.0  mWcm  ym  sr  ,  and  RH=90%,  WS=0,  one  obtains 

% 

HI—  3.8  km.  Again  this  may  be  interpreted  as  a  measure  of 
number  density  in  the  fine  mode  for  a  given  scale  height. 
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8.  SUMMARY  OF  RESULTS 


This  report  has  described  the  technical  approach,  devel¬ 
opment,  and  implementation  of  a  model  to  relate  simulated 
satellite  incident  radiances  (in  the  wavelength  region  from 
0.4  to  1.1  ym)  to  variations  in  meteorological  variables 
such  as  relative  humidity,  windspeed,  aerosol  number  density, 
visual  range,  and  aerosol  scale  height.  In  this  context,  the 
results  of  this  effort  may  be  considered  an  extension  of 
existing  optical  propagation  models  capable  of  relating 
remotely-sensed  and  insitu  optical  transmission  data  bases. 

Technical  elements  treated  in  this  investigation 
included:  (1)  calculation  of  non-aerosol  atmospheric 

transmission  properties  (Rayleigh  scattering  and  gaseous 
absorption) ,  (2)  implementation  of  physical  aerosol  models 

describing  dependence  of  size  distribution  and  index  of 
refraction  on  aerosol  growth  and  wind  speed 
dependence  of  aerosol  oceanic  components,  (3)  calculations 
of  maritime  aerosol  optical  properties  using  Mie  theory  and, 
(4)  incorporation  of  the  above  in  an  efficient  radiative 
transfer  calculation.  The  model  was  exercised  to  provide 
sensitivity  analyses  of  both  wavelength-dependent  and  DMSP 
bandpass  weighted  radiance  to  the  meteorological  variables 


cited  above. 


9.  CONCLUSIONS  AND  RECOMMENDATIONS 


9.1  Conclusions 

A  number  of  conclusions  emerge  based  on  the  results 
reported  in  Section  7. 

.  Calculations  using  both  the  maritime  and  hybrid 
aerosol  models  indicate  that  simulated  emergent 
radiances  in  the  near  infrared  window  regions 
(A>0.7  ym)  are  highly  sensitive  to  the  presence  of 
coarse  mode  (oceanic  origin)  particulates. 

.  Enhanced  radiances  in  these  regions  are  indicated 
for  increasing  relative  humidity  and  windspeed, 
primarily  due  to  the  characteristically  greater 
aerosol  growth  for  these  large  particles  and  the 
mechanical  addition  of  particles  due  to  wind  mixing. 

.  Results  from  the  hybrid  model  indicate  that  the  near 
infrared  effect  is  enhanced  when  higher  wind  speeds 
and  relative  humidity  occur  concurrently. 

.  The  dependence  of  simulated  radiances  in  the  visible 
spectral  region  (A<0.7  ym)  on  these  factors,  by 
contrast,  is  significantly  less  than  that  evidenced 
at  larger  wavelengths.  At  these  shorter  wavelengths 
radiances  appear  to  be  determined  primarily  by 
Rayleigh  scattering  and  accumulation  mode  aerosols. 

.  Calculations  suggest  that  the  DMSP  sensor  is  capable 
of  delineating  reduced  visibility  situations  due  to 


the  dependence  of  visual  range  on  the  number  density 
of  accumulation  mode  aerosols.  However,  a  quantific¬ 
ation  of  the  relationship  between  DMSP  radiances  and 
surface  visual  range  requires  independent  knowledge 
of  other  variables  such  as  relative  humidity  and  fine 
mode  aerosol  scale  height. 

For  fixed  relative  humidities  and  wind  speeds  (if 
applicable) ,  wavelength-dependent,  simulated  radiances 
generally  increase  with  aerosol  scale  height. 

Variation  with  relative  humidity  and  wind  speed  are 
enhanced  with  increases  in  aerosol  scale  height. 

The  effect  of  water  vapor  on  remote  sensing  within 
this  spectral  range  cannot  be  ignored  outside  the 
water  vapor  absorption  bands  due  to  its  coupling 
through  aerosol  growth  laws  with  aerosol  continuum 
scattering  elsewhere. 

Variations  in  ocean  surface  albedo  appear  to  affect 
only  simulated  visible  wavelength  radiances. 

The  dependence  of  simulated  DMSP  bandpass-weighted 
radiances  on  the  meteorological  variables  discussed 
above  suggests  that  the  sensor  may  be  used  to  provide 
correlating  evidence  to  characterize  the  meteorological 
situation  within  the  field  of  view,  given  that  some  of 
the  relevant  variables  are  available  independently  from 
surface  observations  or  other  satellite  sensors. 

Standing  alone,  the  DMSP  sensor  cannot  uniquely 
characterize  the  meteorological  variables  treated  above. 
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9.2  Recommendations 

The  results  of  this  study  may  potentially  provide  an 
economical  tool  to  assess  the  effects  of  meteorological 
variables  on  remotely  sensed  radiances  in  the  visible  and 
near  infrared  portion  of  the  spectrum.  A  broader  applica¬ 
tion  than  the  present  study  concerning  the  DMSP  sensor  may 
be  envisioned  by  examining  other  spectral  regions  and 
incorporating  other  relevant  bandpass  functions.  These 
ambitious  goals  should  be  undertaken  with  great  caution, 
however,  until  a  number  of  matters  arising  from  the  present 
study  are  resolved. 

A  variety  of  recommendations  can  be  made  concerning  the 
physical  elements  of  the  modeling  approach  adopted.  These 
include:  (1)  examination  of  other  aerosol  growth  laws  and 
physical  models,  (2)  implementation  of  a  surface  albedo 
model  coupling  surface  windspeed  to  sun  glitter  effects,  and 
(3)  improved  treatments  of  radiative  transfer.  As  a  minimal 
requirement,  with  respect  to  the  last  item,  results  reported 
here  should  be  compared  with  more  accurate,  numerical 
radiative  transfer  treatments  run  with  the  same  data  set  of 
optical  parameters  to  better  quantify  errors  associated  with 
the  approximate  method  adopted. 

This  study  was  conceived  and  initiated  as  a  theoretical, 
modeling  exercise.  The  utility  of  the  approach  presented  and 
the  validity  of  the  results  reported  must  be  verified  by 
conducting  a  validation  study  based  on  an  appropriate  set 
of  field  data  including  aerosol  optical  and  physical 
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parameters,  meteorological  variables,  and,  if  possible, 
concurrent  remotely-sensed  radiances. 
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12.  APPENDIX:  DERIVATION  OF  APPROXIMATE  ANALYTICAL 

SOLUTION  TO  RADIATIVE  TRANSFER  EQUATION 

Following  Kaufman  (1979) ,  the  analytical  solution  to  the 
radiative  transfer  Equation  [  (33)  —  (34) ]  used  in  this  study  is 
derived  by  approximating  the  integral  term  in  the  source 
function  J  [Equation  (34) ]  by  assuming  that  the  upward  (y  >  o) 
and  downward  (y  <  o)  diffuse  (i.e.  scattered)  intensity 
fields  are  independent  of  direction  (y,<J>)  and  given  by  I+ 
and  I~,  respectively.  Mathematically,  this  approach  may  be 
considered  as  a  two-point  quadrature  of  the  angular  integra¬ 
tion  over  upper  and  lower  hemispheres.  Physically,  the 
method  is  interpretable  as  an  approximation  of  the  multiply 
scattered  contributions  to  the  intensity  field.  For  incident 
solar  irradiance  ttF,  Equation  (34)  is  approximated  as: 

Fe-T/yo 

Js(x,+  y,<j>)  m  j-  P(y,+  y0?  4>f<t>0) 

1  2tt 

+  1  1  P(y,±tf;4>r4>  ')  dy'd4>' 

o  o 

1  2  IT 

+  |  P(]i/  +  y';  dy  *d<j) '  (A .  1 ) 

0  0 

Upon  substitution  of  (A.l),  the  radiative  transfer  Equation 
(33)  is  written  separately  for  each  of  the  two  streams,  I—, 
as : 
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2tt  1 


dl 


=  I+  -  j  J  P(y,  +  y' ;  <t>,<t>')  dy'd<J>' 


o  o 


2  IT  1 

£  /  /  P(M 


-  y  '  ;  (J)  ,4) ' )  dy  'd<j> ' 


o  o 


F  t,  ,  ,  ,  .  -T/yo 

j  P  (y  ,  -  y0  ;  <J>  ,4»0  )  e 


(A. 2a) 


2tt  1 


dl 

y  3 —  =  I 
dr 


~  T¥  I  J  P(u'  "  y,;  <*>'<l>,)  dy’  dct> » 


0  o 


_  2tt  1 

~  hr  J  J  p(p'  y,;  ^ ^ ^  dlJ’  d<f)' 


o  o 


^P(y,  y0  ;  <}>,<J>0)  e  T</y° 


(A. 2b) 


The  angular  integrations  in  (A.2a,b)  are  performed  by  noting 
the  phase  function  normalization: 

,  2  ir  +1 

U)°  =  4?f  J  P(y'  y,;  *'  *’)  dy’  d*’  (A-3) 

o  -1 

and  defining  the  azimuthally  symmetric  function,  P°;  the 
backscatter  fraction  for  monodirectional  radiation,  8;  and 
the  backscatter  fraction  for  isotropic  radiation,  $’  as: 
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r 


P° (y  »u ' )  =2 7  J  P(y,y ' ; 4> , <J> ’ )  d$' 
o 

1 

6  (y)  =  2^  /  p°<^'>  dy’ 

o 

1 

B ’  =  J  B(y’)  dy’ 
o 


Applying  these  definitions,  the  intensity  Equations 
become : 

y  =  i+  -  i+  0)O  [l-B(y)]  -  l“w0B(y) 

-  |  P(y r -y o ? <l> r 4>0 )  e  T/y° 

-  y  Jp-  =  I-  -  I+  0)0B  (y)  -  I"  o)0  [  1-S  (y )  ] 

-  ^  P(y ,y0;4>f<f>0)  e  T//y° 


Assuming  that  the  scattered  diffuse  intensities  are 
by  the  fluxes  F^t)  such  that: 

I*  (x  ,y  ,<J>)  =  F±(t)/u 


(A. 4) 


(A. 5) 


(A. 6) 


(A. 2) 


(A. 7a) 


(A. 7b) 


given 


(A. 8) 
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Equation  (A. 7)  may  be  written  as: 


y 


TT 


F+[l-8(y) ]  +  F“B(u) 


-  ^  P(y#-yo«<J>»4>o)  e  T^y° 


(A. 9a) 


-  y 


dl 

dx 


I 


u>o 

IT 


F+3  (y)  +  F_  [ 1-6 (y ) H 


-  P  (y /y0 ; <)> »<j>0)  e 


-T/yc 


(A. 9b) 


In  matrix  rotation,  these  correspond  to  equations  (40),  (41), 
and  (42)  in  the  text. 

The  flux  equations  (equations  (35),  (36)  and  (37)]  in 
the  text  may  be  obtained  by  integration  of  equations  (A. 9) 
over  the  angle  y  to  give: 


dF+ 

37"  " 


2F 


U-Wod-B')  ] 


2F  w0B ’ 


ttFw06  (y0) 


e-T/wo 


(A. 10a) 


dF 

dx 


-2F-  [  l-ti)0  (1-B  * )  ]  +  2F+  u>0(5  ' 


tiFw0(  1-6  (y0)  ] 


e"T/P(> 


(A. 10b) 


Applying  the  boundary  conditions  [of  equation  (38)  in  the 
text] : 
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r 


F  (t=o)  =  0 


F+(t=t*)  =  Aj^  [7T  y0  e 


-T  /y°  +  F~ (  T  =  T  * )  ] 


(38) 


the  solutions  for  F—  are  [cf  (39)  in  text] : 


F+(t)  =  AekX  +  Be“kx  +  Ce-T^Wo 


F  (T)  = 


^2 


A(Yi"k)ekx  +  B(Y1+k)e-kT+  Ye-x/,y° 


(A. 11a) 


>  (A. lib) 


where  the  appropriate  constants  are  given  by: 


y1  =  2  [l-u0(l-B')] 


(A. 12) 


Yo  =  2  (1)08* 


(A. 13) 


k  =  (Yj2  *  Y22)'i 


(A. 14) 


A  =  [BfYj+k)  +  Y]/(k-Y1) 


(A. 15) 


B 


=  (E-^1”  +  E2e“T  ^yo)/(E3ekx  +  E4e_kx  )  (A. 16) 


C  =  TT  U)0 


-  y  I)  ( y 0 )  -  y2[1-6(Vo>  ] 

Uo  1 


■(ot) 


(A. 17) 
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(A. 18) 


f 


Y  =  C  (Yi  +  )  -  ttF  u)06(y0) 


Ex  ■  Y  [l/(Y1"k)  -  A2/y2 1 

(A. 19) 

aly 

E2  =  [-  c  +  nFy0  Al  +  -y-  ] 

(A. 20) 

E3  =  (Yj.  +  k)[l/(k-Yl>  +  AL/Y2] 

(A. 21) 

E4  =  [1  -  Al  (Yi+k)/Y2l 

(A. 22) 

The  solution  for  the  upward  (i.e.  emergent)  intensity, 
I+  (t  ,y  ,<J>)  ,  is  obtained  by  integrating  equation  (A.  9)  subject 
to  the  boundary  conditions  [see  text  (43),  and  (44)]: 


I  (t=o,  y,<J>)  =  0 


(43) 


I+  (T=T*,M,<j>) 


F+  (t=t*)/ti 


(44) 


to  obtain  [compare  to  text  (45)]: 


0,ekT  U,e'k1  u3e-T/l,o 


I-  * . .  -  Dl  +  ♦  ^7^  «»•»> 


where : 


Di  = 


F*  (  T*  )  _ 
IT 


UiekT  U2e"kT 


1-yk  1+yk 


_  U3e 


-x*/y0 


i+y/y< 


-T*/y 


(A. 24) 
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UL  =  HT  tl-6(u)  +  (Yi-k)]  (A. 25) 

u2  =  ~  U-0(M)  +  (Yl+k)J  (A. 26) 

Cu>0  p  . 

U3  =  —  [1-6  (u)  ]  +  —  SIHL  Y  +  |  P  (|I,-P0;*,*0)  (A. 27) 

The  required  value  of  F  (x  )  in  D3  (A.24)  is  available  from 
equation  (A.  11a).  The  values  of  B',3(v0»  and  3(h0)  in  the 
above  expressions  are  obtained  from  Figure  22,  after  calcu¬ 
lating  the  appropriate  asymmetry  parameter,  g,  using 
equation  (48)  .  Values  of  x*  and  u)0  used  above  are  obtained 
from  equations  (46)  and  (47)  ,  respectively,  using  the  modeled 
optical  data  for  gases  and  aerosols  described  in  Sections  3 
and  5,  respectively. 
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